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Abstract: Precision studies of the production of a high-transverse momentum lepton 
in association with missing energy at hadron colliders require that electroweak and QCD 
higher-order contributions are simultaneously taken into account in theoretical predictions 
and data analysis. Here we present a detailed phenomenological study of the impact of 
electroweak and strong contributions, as well as of their combination, to all the observables 
relevant for the various facets of the pp — > lepton + X physics programme at hadron 
colliders, including luminosity monitoring and Parton Distribution Functions constraint, 
W precision physics and search for new physics signals. We provide a theoretical recipe to 
carefully combine electroweak and strong corrections, that are mandatory in view of the 
challenging experimental accuracy already reached at the Fermilab Tevatron and aimed at 
the CERN LHC, and discuss the uncertainty inherent the combination. We conclude that 
the theoretical accuracy of our calculation can be conservatively estimated to be about 2% 
for standard event selections at the Tevatron and the LHC, and about 5% in the very high 
W transverse mass/lepton transverse momentum tails. We also provide arguments for a 
more aggressive error estimate (about 1% and 3%, respectively) and conclude that in order 
to attain a one per cent accuracy: 1) exact mixed 0{aa s ) corrections should be computed 
in addition to the already available NNLO QCD contributions and two-loop electroweak 
Sudakov logarithms; 2) QCD and electroweak corrections should be coherently included 
into a single event generator. 
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1. Introduction 

The electroweak gauge bosons and Z° were discovered at the SpS collider at CERN 
more than twenty years ago Q. Albeit a long time has passed since then, the production 
of electroweak gauge bosons in hadronic collisions is still a topic of deep interest in mod- 
ern particle physics. Actually, single W and Z boson production is used at the Fermilab 
Tevatron collider to derive precision measurements of the W boson mass and width (2[, 
to extract the electroweak (EW) mixing angle from the forward-backward asymmetry in 
the neutral current (NC) channel || and to severely constrain the Parton Distribution 
Functions (PDF) through the measurement of the W charge asymmetry These pro- 
cesses are also used for detector calibration and the measurement of their total production 
cross section can be compared with the corresponding QCD prediction, in order to test the 
convergence of the strong coupling expansion in perturbative QCD calculations ||. 
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In the near future, at the Large Hadron Collider (LHC) at CERN, the production of 
W and Z bosons will continue to be a relevant process, because of its large cross section 
and very clean signature given by one isolated charged lepton with missing transverse 
energy (for W production) and by two isolated charged leptons with opposite charges (for 
Z production). In particular, thanks to the very large statistics, a measurement of the W 
mass with an uncertainty of about 10 MeV should be feasible at the LHC [||, ||. These 
processes are also good candidates to understand the detectors performances in the early 
stage of data analysis at the LHC JnJ, to monitor the collider luminosity with per cent 
precision |ll], 12, and constitute a background to new physics searches, noticeably new 



heavy gauge bosons Z 1 and W [[H], 14], whose discovery is an important goal of the LHC 

Hi- 

It is important to realize that the many facets of the W/Z physics programme at hadron 
colliders require, for obvious reasons, the measurement of different observables, depending 
on the physics goal of interest. Correspondingly, a number of observables must be precisely 
predicted and simulated, to avoid theoretical bias in data analysis. For example, for precise 
W mass measurements the relevant observables are the W transverse mass and the lepton 
transverse momentum, while luminosity studies require a deep understanding of the total 
cross section, the W/Z rapidity and lepton pseudorapidity. On the other hand, the lepton 
pair invariant mass produced in Z production and the W transverse mass in their high 
tails are the observables to be focused on for the search for new physics signals. 

To fully exploit the potential of the Tevatron and the LHC for all the above physics 
goals, the theoretical predictions have to be of the highest standard as possible. In partic- 
ular, the high luminosity at the LHC implies that systematic errors, including theoretical 
ones, will play a dominant role in determining the accuracy of the cross sections. This 
requires to make available calculations of W and Z production cross sections including 
higher-order corrections originating from the EW and QCD sector of the Standard Model 
(SM). Furthermore, the implementation of such calculations in Monte Carlo (MC) gen- 
erators is mandatory, in order to perform realistic studies of the impact of higher-order 
corrections on the observables and to compare theory with data. 

In this paper we study the impact of electroweak and strong corrections, as well as of 
their combination, to the production process pp — > W — > lepton + X. We also include 
in our analysis photon-induced contributions to W — > lepton + X production and we 
provide results for top-pair production process, which may give rise to the signature of a 
high-transverse momentum lepton in association with missing energy, in order to include 
a non-negligible contribution for a precise evaluation of the overall background to new- 
physics searches. A similar study, which some of us have coauthored, has been recently 
performed for the NC Drell-Yan (DY) channel in the high invariant mass tail in p5| , 
while an evaluation of the theoretical uncertainties in the cross sections of the charged 
current (CC) and NC processes around the W/Z peaks has been addressed in []16|, |l7j by 
strictly following the procedure presented by our group in some conference proceedings [18]. 



However, it is important to notice that at a difference with respect to (w.r.t.) the works [16 



17], where the effects of EW and QCD corrections are separately studied, in the present 



paper the combination of such effects is analyzed in detail, by providing different theoretical 
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recipes and discussing the related uncertainties. Furthermore, we provide also predictions 
for the transverse mass region above 1 TeV, which is an important range for new physics 
searches. 



As far as complete O(a) SM EW corrections to pp — > W — >• lv\ are concerned, they 
have been computed independently by various authors in ^0], 22, 23, 24, [25|. Also 
one-loop corrections in the Minimal Supersymmetric Standard Model have been calculated 
|2l|. The predictions of a subset of such calculations have been compared, with the 



m 



same input parameters and cuts, in [27, 28], finding a very satisfactory agreement. The 
contributions due to the emission of multiple photon radiation for the observables of interest 
for W precision physics have been calculated in |2^. 



31 , p2| , |33| , Eg]. Comparisons of 



such multi-photon calculations can be found in [32, 33, p4l , showing good agreement. 



Concerning QCD calculations and tools, the present situation reveals quite a rich 
structure, that includes NLO/NNLO corrections to W total production rate pBL |3q| , fully 
differential cross sections known at NLO/NNLO QCD |37|, ||, ||, go], gl], [42|, resummation 
of leading/next-to- leading logarithms arising from soft gluon radiation flij ], NLO cor- 



rections merged with parton shower evolution [45, Bq, 47], as well as leading-order matrix 



while 



elements generators matched with vetoed parton shower [49|, |(], 51]. 

A review of the theory of W and Z boson production can be found in |5^, [_ 
recent developments about QCD and generators are described in |54|, |55| . 

On the other hand, it must be noticed that the combination of EW and QCD cor- 
rections is still at a very preliminary stage, only a few attempts being known in the liter- 
ature [^, 57|. In particular, in |)6| ] the effects of QCD resummation are combined with 
the leading part (in the statistically dominant region of the process around the W mass 
peak) of EW corrections, given by NLO final-state QED radiation and improved Born 
approximation to account for the pure weak corrections. Therefore, in the latter analysis 
the effect of purely (Sudakov-like) EW logarithms is missing, when it is known that the 
full set of EW corrections is required in addition to photonic corrections to reach a per 

the treatment of QCD effects 



cent theoretical precision 21, [2^, |2q] . Furthermore, in 
is limited to the inclusion of multiple soft-gluon corrections, leaving room to more detailed 
studies of the QCD sources of contributions. The present situation of EW calculations 
and the recent development of more sophisticated QCD generators makes feasible a more 
complete analysis, both for what concerns QCD and EW effects. 
Therefore, the aims of the present paper are the following: 



1. to investigate the various sources of QCD effects (PDF uncertainties, comparison 
of presently available perturbative QCD calculations, allowing for variations of the 

to consider 



factorization/renormalization scale), extending the analyses of [40, 56 



the full set of observables of experimental interest and to make use of computational 
tools which include recent advances in QCD phenomenology; 

to combine the results of the state-of-the-art of EW calculations with the QCD cal- 
culations and generators, that are typically used at the Tevatron and the LHC, in 
order to investigate the role of EW effects in association with QCD contributions; 
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3. to assess the reliability of existing tools and demonstrate the need of a deep un- 
derstanding of QCD and of EW-QCD combination for precision physics studies at 
hadron colliders. 



The outline of the paper is as follows. In Section ^ we review the main features of the 
theoretical tools used in our study, describing the treatment of PDF uncertainties (Section 
|2.1| ), the implementation of QCD effects according to different degrees of approximation 
in existing generators (Section |2.2|) , the treatment of NLO EW corrections and photon- 
induced processes (Section |2.3| ), as well as the combination of EW effects with strong 
corrections (Section |2.4[ ), which is the main theoretical feature of the paper. In Section 
|H we present and discuss the numerical results for the observables that are grouped on 
the basis of their specific experimental motivation, both at the energies of the Fermilab 
Tevatron proton-antiproton (pp) collider and of the CERN LHC proton-proton (pp) collider. 
Conclusions and possible developments of our work are drawn in Section ||. 



2. Theoretical ingredients 

2.1 Parton Distribution Functions (PDF) 



It is widely documented [59| that a precise knowledge of the partonic structure of the pro- 
ton will be an essential ingredient for the physics potentials of the LHC, both for what 
concerns discovery and precision physics. At present, there is a great deal of interest in 
understanding how the uncertainties on the determination of the Parton Distribution Func- 
tions (PDF), both of experimental and theoretical origin, translate into uncertainties for 
those observables that will be measured at the LHC. A short review of recent progress 
in Parton Distributions and implications for LHC physics can be found in |30|. A par- 
ton distribution library, known as LHAPDF is available and gives the possibility of 
comparing different PDF parameterizations. 

In particular, some sets allow to estimate how the errors of the experimental measure- 
ments affect the PDF parameterization within a certain practical tolerance. The latter is 
defined as the maximum allowed of the Ax 2 variation w.r.t. the parameters of the best 
PDF fit. In our analysis, among the PDFs sets available in LHAPDF, we will make use of 
MRST2001E |6^] and CTEQ6 Q, as done in previous similar studies [64|]. For these 
two different sets of NLO PDFs the tolerance T is assumed to be T = (MRST2001E) 
and T = 10 (CTEQ6), resulting in PDFs uncertainties larger for a factor of v2 for CTEQ6. 
It is worth emphasizing that the uncertainties obtained according to such a procedure are 
purely experimental (i.e. as due to the systematic and statistical errors of the data used 
in the global fit), leaving aside other sources of uncertainty of theoretical origin, such as, 
for example, uncertainties due to the truncation of DGLAP perturbation expansion, choice 
of the parameterization of non-perturbative input PDF or ln(l/x) effects. However, they 
can be estimated separately and it is known from previous studies |35|, ^36| that the theory- 
induced uncertainty is similar in magnitude to the uncertainty due to experimental errors 
at the Tevatron, but larger at the LHC. However, in analogy to previous phenomenological 
studies p8|, |64[|, we limit to consider PDFs uncertainties deriving from errors on the data, 
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just to get an idea of the order of magnitude of such effects in comparison with the size of 
perturbative corrections that are the main subject of our analysis. 

Finally, we point out that, since we include in our predictions EW corrections and, 
therefore, we need to treat collinear singularities due to photon radiation off the initial 
state quarks for the sake of consistency at NLO EW, it is necessary a set of PDF including 
QED effects in DGLAP evolution. The MRST group performed a global PDFs analysis 
including QED contributions, which have been incorporated in the set of PDFs known as 



MRST2004QED [67]. It turns out that the QED evolution slightly modifies the standard 



QCD-evolved PDFs for typical x and Q 2 values probed by weak boson production at the 



Tevatron and the LHC, as already argued in previous studies |6^, 70]. However, it 
must be stressed that the inclusion of QED corrections into PDFs dynamically generates 
a photon density inside the (anti)proton. This implies that EW processes initiated by a 
photon (photon-induced processes) contribute to EW gauge boson production in hadronic 



collisions, as discussed in Section |2.3| . As largely discussed in the literature, the present 
PDF uncertainties are expected to be significantly reduced when LHC data will become 
available. 

The parameterizations MSTW [JFl] and CTEQ6.6 |72||, appeared in the literature dur- 



ing the completion of the present work, will be discussed shortly in Section 3.4. 



2.2 QCD corrections and Monte Carlo tools 

The QCD codes typically used in high-energy collider experiments are "general purpose" 
parton shower MC's, such as, for example, HERWIG @ or PYTHIA Q, Apacic++ |7§ 
or MC programs based on a fixed-order perturbative calculation of a given process. How- 
ever, the new challenges of the Tevatron Run II and the LHC have spurred on improved 
approaches to QCD MC's and calculations, including the matching of exact leading-order 
matrix elements for multiparticle production with parton shower, the matching of NLO 
corrections with parton shower, advances in the techniques for NNLO calculations, as well 



as improvements in the parton shower algorithms [54, 77, 78 1. 

As already remarked in the Introduction, we make use in our analysis, for the sake 
of definiteness, of those QCD calculations and relative codes that are already used at the 
Tevatron or will be presumably used at the LHC for data analysis, and, at least in our 
opinion, are representative of those advances sketched above. We therefore consider the 
following set of programs 



1. ALPGEN |4S]: it is a code for the generation of many multiparton processes in 
hadronic collisions, computing exact matrix elements by means of the ALPHA algo- 
rithm [79]. The leading-order matrix elements are matched with appropriate Sudakov 



form factors to add parton shower effects (according to the so-called MLM prescrip- 



tion to match matrix elements and parton shower [pQ , pl| , |82 , 83| , |84J ) and by vetoing 
shower evolution leading to multiparton final states already described by the matrix 
element computation. A review of the matching algorithm implemented in ALPGEN 
and here used can be found in [ 84 1 , while a comparison between the predictions of the 
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different procedures of matrix elements-parton shower matching is given in |85|, [Sq] . 
Additional information is available at http://mlm.home.cern.ch/mlm/alpgen/ 

2. MCFM [ffi?f| : it is a program for the MC simulations of various processes at hadron 
colliders, giving predictions, for most processes (including W/Z production), at NLO 
in QCD and incorporating full spin correlations. It is available, together with a 
detailed documentation, at http://mcfm.fnal.gov 



3. MC@NLO ||45|| : it is a generator incorporating NLO QCD corrections into the parton- 
shower generator HERWIG for a large class of processes. As such, MC@NLO pro- 
vides a consistent description of exact 0(a s ) emission effects together with a leading- 
logarithmic resummation of soft and collinear QCD radiation. More details can be 
found at http://www.hep.phy.cam.ac.uk/theory/webber/MCatNLO/ 



4. ResBos [43 , 44]: it is a MC integrator program, presently used at the Tevatron, that 
computes fully differential cross sections of various boson production processes in 
hadron-hadron collisions, including DY-like processes, either with NLO initial-state 
QCD corrections, or with soft-gluon initial-state resummed QCD corrections, the 
latter according to the Collins, Soper and Sterman (CSS) resummation formalism 
Recently, Resbos has been improved by the inclusion of final-state NLO electro- 
magnetic contributions to W boson production [56], resulting in the ResBos-A code. 



More information can be found at http://hep.pa.msu.edu/resum/. Since in our 
study we are interested in both QCD and EW effects, we use, among the different 
ResBos versions available on the web, the ResBos-A generator, including CSS resum- 
mation and final-state QED radiation, making use of the grids publicly available for 
the Tevatron energy 



5. FEWZ [40, |4l[: it is a code computing fully differential cross sections at NLO/NNLO 



in perturbative QCD for W and Z boson production in hadron collisions, including 
finite-width effects and full spin correlations. 

It is available at http://www.phys.hawaii.edu/~kirill/FEHiP.htm 

As far as FEWZ is concerned, it will be used in the present paper for tuning purposes at 
the NLO level for fully integrated cross sections. Its predictions at the NNLO level and the 



phenomenological relevance of NNLO corrections have already been analyzed in (T^, 17]. 

It should also be mentioned that a recent advance in the QCD predictions for DY 
processes is represented by the calculation of f46| . This consistently combines the NLO 
QCD corrections with parton shower resummation, avoiding negative weights in the gen- 
eration of events, and is available in the MC POWHEG. We will not show in our analysis 



predictions of POWHEG, since they are already compared with MC@NLO 47], finding 
fair agreement for j'ei-inclusive quantities. 

2.3 Electroweak contributions 

The radiative corrections originating from the EW sector of the SM are taken into ac- 
count in our analysis according to the recent calculation of [25], which the reader is ad- 
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dressed to for more details. At a variance of QCD corrections, we do not need to take into 
consideration here the results of different, independent calculations because it is known 
from recent tuned-comparison studies p7| , f2§| 1 that the predictions of those calculations 
typically used or under consideration by hadron collider collaborations, such as DK pl|| , 
SANC and WGRAD [^C], are in very good agreement with the results of HO- 
RACE, which is the generator implementing the theoretical formulation of available 
at http://www.pv.infn.it/~hepcomplex/horace.html. In HORACE, the exact one- 
loop EW corrections to the CC DY process are consistently matched with higher-order 
leading logarithms due to multiple photon emission, that have been shown to be not ir- 
relevant for the precision target of W mass measurements at hadron colliders [29, BfJ. By 
itself, HORACE does not include any effect due initial-state QCD parton shower, but the 
events generated with it can be interfaced with standard shower MC programs. 

Since we are interested in a detailed description of EW effects too, we need to include 
in our calculation, for the sake of consistency, also the contribution of the so-called photon- 
induced processes -yq — > q'lvi and jq 1 — ► qlvi (see Figure ||). Actually, these processes are 
of the same EW perturbative order of one-loop EW corrections and give rise to the same 
signature of the inclusive W production process. Furthermore, their contribution has been 




w+ 

d X i + d 

Figure 1: Examples of Feynman diagrams for the photon-induced processes 




shown by Dittmaier-Kramer in [27] to be non negligible for VF-boson production in the 
presence of realistic event selection cuts at the LHC. The matrix elements for the photon- 
induced processes (73 — > q'lvi) are obtained by crossing symmetry from the qq' — > Iviy 
matrix elements, already available from the calculation of hard-bremsstrahlung diagrams 
contributing to EW corrections. 

It is worth noting that the inclusion of photon-induced processes introduces in the 
0(a)-corrected parton-level cross section mass singularities due to a collinear splitting 
7 — * QQ- As for photon radiation off the initial-state quarks present in the NLO EW 
calculation and, again, in analogy with factorization in NLO QCD calculations, the collinear 
singularities due to 7 — > qq collinear splitting can be absorbed in the distribution functions, 
by replacing the (anti-)quark distribution as described in detail, with all the necessary 
formulae, in flSST . 
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2.4 Combination of electroweak and QCD corrections 

The aim of this Section is to describe how the contributions of EW origin are combined, 
in our approach, with QCD predictions. 

A first strategy for the combination of EW and QCD corrections consists in the fol- 
lowing additive formula 



da 


J da \ J 


da 




da 




dO 


QCD&EW I dO J q CD [ 


dO_ 


EW 


dO 


loIherwig PS 



(2.1) 



where da/dOQ CD stands for the prediction of the observable da/dO as obtained by means 
of one the QCD generators described in Section [2]^, da/dO EW is the HORACE prediction 
for the EW corrected da/dO observable, and da/dO LO is the lowest-order hadron-level 
result. The label HERWIG PS in the second term in r.h.s. of eq. ( [2.1D means that EW 
corrections are convoluted with QCD PS evolution through the HERWIG event generator, 
in order to (approximately) include mixed 0(aa s ) corrections and to obtain a more realistic 
description of the observables under study. In principle, concerning QCD, the predictions 



of any of the codes considered in Section 2.2 could be used in eq. ( |2.1| ). However, because in 
our approach EW corrections are treated at NLO and are convoluted with the QCD parton 
shower cascade by HERWIG, it follows that, in eq. (2.1), the QCD contributions can be 
consistently combined with EW effects only by making use of the results provided, among 
the different QCD generators, by MC@NLO, in order to treat QCD NLO and shower cor- 
rections on the same ground of the EW counterpart. Therefore, in the following, we will 
show results for the QCD and EW combination assuming [da / dO]Q CD = [(i<r/(iC'] MC@NL0 1 . 
This implies, in turn, that the HORACE and MC@NLO generators need to be properly 
tuned at the level of input parameters and PDF set, in order to provide meaningful predic- 
tions. Furthermore, since we are interested in consistent comparisons between the results 
of independent codes, this tuning procedure has to be applied to each generator in use, as 
discussed in Section 3.2. Concerning the convolution of NLO EW corrections with QCD 
PS, it is worth noting that, according to eq. (2.1), the contributions of the order of aa s 
are not reliable when hard non-collinear QCD radiation turns out to be relevant, e.g. for 
the lepton and vector boson transverse momentum distributions in the absence of severe 
cuts able to exclude resonant W production. In this case, a full 0(aa s ) calculation would 
be needed for a sound evaluation of mixed EW and QCD corrections. Full 0(a) EW cor- 
rections to the exclusive process pp — > W + j (where j stands for jet) have been recently 



computed, in the approximation of real W bosons, in |p9| , 90 1, and for off-shell W in fl91fl , 
while one-loop weak corrections to Z hadro-production have been computed, for on-shell 
Z bosons, in p2| . It is also worth stressing that in eq. ( |2.1| ) the infrared part of QCD 
corrections is factorized, whereas the infrared-safe matrix element residue is included in an 
additive form. 



1 Alternatively one could use POWHEG in association with PYTHIA PS, without spoiling the general 
procedure presented here. 
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It is otherwise possible to implement a fully factorized combination as follows: 



_ | -j, [dv/dO] MC@NLO - [da / dO] RERW1G PS \ ^ 
qcd®ew \ [d a /dO] LO / NLO J 

L «^EW J HERWIG PS 



da 
dO 



where the ingredients are the same as in eq. fl2.1| ) but also the QCD matrix element 
residue in now factorized. It is worth noticing that the QCD correction factor in front 
of {da/dOEw} herwig PS ^ s defined in terms of two different normalization cross sections, 
namely the LO or the NLO one, respectively. As it can be easily checked, the two prescrip- 
tions differ at the order a 2 by non-leading contributions. Nevertheless, eq. ( |2.2[ ) normalized 
in terms of the LO cross section can give rise to pathologically large order a 2 corrections 
in the presence of huge NLO effects, as will be discussed in the next section. On the other 
hand, when NLO matrix element effects do not introduce particularly relevant corrections, 
the two prescriptions are substantially equivalent, as it will be demonstrated too in the 
discussion of numerical results. Some comments are in order here. Equations ( [2.1D and 
(|2.2| ) in both prescriptions have the very same 0(a) and 0(a s ) content. While the two fac- 
torized prescriptions contain the same 0(aa s ) contributions, the additive prescription has 
mixed 0(aa s ) contributions which differ from the factorized recipes. Therefore, their rela- 
tive difference can be taken as an estimate of the uncertainty of QCD and EW combination 
and will be discussed in detail in the next section. 

3. Numerical results and discussion 

3.1 Input parameters and event selection criteria 

The numerical results have been obtained using in each code the following values for the 
input parameters 





1.16639 10" 5 GeV~ 2 


m w 


= 80.419 GeV 


m z 


= 91.188 GeV 


T w = 


2.048 GeV 


■ 9 

sin 


9 w = l- m 2 w /m 2 z 


m H 


= 120 GeV 


m e = 


510.99892 KeV 


"V 


= 105.658369 MeV 


m T 


= 1.77699 GeV 


m u = 


320 MeV 


m c 


= 1.55 GeV 


m t 


= 174.3 GeV 


m d = 


320 MeV 


m s 


= 500 MeV 


m b 


= 4.7 GeV 


v cd = 


0.224 


Vcs 


= 0.975 


Vcb 


= 


v ud = 


V 1 - V cd 


Vus 


= 0.224 


Vub 


= 


Vtd = 





V ts 


= 


Vtb 


= 1 



and adopting the input scheme for the calculation of EW corrections, where, in particu- 
lar, the effective electromagnetic coupling constant is given in the tree-level approximation 
by 



tree _ Vjjgjf sin2 
a G u ~ - lo-J-J 
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However, for the coupling of external photons to charged particles needed for the evaluation 
of photonic corrections we use a(0) = 1/137.03599911. 

We study the production process pp — * lepton + X, both at the Tevatron pp collider 
(y/s = 1.96 TeV) and at the pp LHC (y/s = 14 TeV), presenting results only for final-state 
muons, for the sake of definiteness. Actually, the difference between electron and muon 
final states is confined at the level of pure EW and multiple photon corrections, and can 
be inferred from the literature addressing such a topic in detail pO, 21, 25, 2S, 30]. 



Tevatron 



p\ > 25 GeV $ T > 25 GeV and \rji\ < 1.2 
pf < 50 GeV 50 GeV < Afj„ < 200 GeV 



Table 1: Selection criteria for the Tevatron 



LHC 

a. p l ± > 25 GeV $ T > 25 GeV and \rji\ < 2.5 

b. the cuts as above © Mf > 1 TeV 

Table 2: Selection criteria for the LHC 



The selection criteria used for our analysis at the Tevatron and the LHC are summa- 
rized in Table || and Table ||[ respectively. In the Tables p 1 ^ and rji are the transverse mo- 
mentum and the pseudorapidity of the muon, Ifi? is the missing transverse energy, which we 
identify with the transverse momentum of the neutrino, as typically done in phenomeno- 
logical studies [2C, 21, 25, 56, 58 1; p^ and are the transverse momentum and 



transverse mass of the W boson, M\ v is the invariant mass of the muon-neutrino pair. The 
selection criteria in the first line of Table |l| are introduced to model the acceptance cuts 
used by the experimental collaborations at the Tevatron, while the constraints on p^ and 
v in the second line are imposed in order to include in our study of the QCD corrections 
the predictions of ResBos-A vl.l, which makes publicly available grids for numerical inte- 
gration corresponding to such constraints. For the LHC, we consider two different set up, 
labeled as a. and b. in Table ^. The first one corresponds to typical cuts used for LHC 
simulations, while for the second one a severe cut on the W transverse mass is superim- 
posed to the cuts of set up a., in order to isolate the region of the high tail of M|f , which 
is interesting for new physics searches. 

In order to avoid edge effects in the final phenomenological analysis, introduced by cuts 
at the generation level, we make use of the following loose generation cuts p l ± > 5 GeV, 
]pr > 5 GeV, \t]i\ < 5 (for both the Tevatron and set up a. at the LHC), with, additionally, 
MY > 0.8 TeV for set up b. at the LHC. These cuts have been imposed simply to improve 
the efficiency in the generation of MC events. For the ALPGEN code, we also use the 
following generation cuts for additional QCD partons: px > 20 GeV, \r]j\ < 5, ARij > 0.7, 
where ARij is the separation between the i-th and j-th parton in the f] — <f> plane. The 
parameters required in ALPGEN for the matching between matrix elements and shower 
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evolution are chosen as: £r,cZus. = 25 GeV, \rjj\ = 5, Rj = 1.05, as done in previous studies 
1 84, 85, |86|. A first attempt to assess the intrinsic systematics associated with the matching 
procedure has been presented in [86]. For our purposes we assume that the up to one jet 
samples are well described by NLO and next-next leading logarithmic codes. Therefore 
the impact of the uncertainties of the matching procedure is related to the importance of 
higher parton multiplicities for the various distributions under consideration. As it will be 
shown in the following, their size is at the level of (or below) 1% at the Tevatron and the 
LHC (event selection a), while they become more relevant at the LHC, event selection b, 
where a more thorough investigation could be worthwhile. 

The set of parton density functions used in our study is CT6EQ6M |63| for the case 
of the Tevatron, with factorization/renormalization scale hr = fip = Mi v for both the LO 
and NLO predictions. Again, this choice has been dictated by the necessity of including 
ResBos-A in our comparisons. Although it is theoretically inconsistent using a PDF set, like 
CT6EQ6M, without QED contributions when combining QCD corrections with NLO EW 
contributions, it is also known, as previously discussed, that the effect of QED evolution 
is small and translates into negligible numerical effects on the observables of interest here 
within the adopted selection criteria. 

For the LHC, all the analysis is performed using the set MRST2004QED, in order 
to consistently incorporate EW corrections and photon-induced processes in association 
with QCD corrections. For the LHC, the QCD factorization/renormalization scale and the 
analogous QED scale (present in MRST2004QED) are chosen to be equal, as usually done 



in the literature |20|, |2l], [22|, |25|, and fixed at (j,r = hf = y P±w + ^lu' wnere Mi u , as done 
in previous studies [58], both for LO and NLO results. 

The numerical results showing the estimate of the PDF uncertainties on the observables 
at the LHC energies (Section |3.4| ) have been obtained using the generator MC@NLO version 
3.3 implementing the LHAPDF package. Results for the NLO MRST2001E (with 30 error 
sets) and CTEQ6 (with 40 error sets) parameterizations as available in LHAPDF are shown, 



using as default scale choice fiR = [ip = yP±\y + 

For a s (Q 2 ) we use as input parameters a s (Mz) = 0.118 at the Tevatron and a s (Mz) = 
0.121 at the LHC. We checked that for all the codes under consideration the values of a s 
agree at the 0.1% level when varying Q from Mz to a few TeV. 

All the numerical simulations have been done without taking into account the effects 
of the underlying event and hadronization processes. 

The numerical results discussed in the following have been obtained according to a 
MC statistics of approximately 10 6 - 10 8 events for the integrated sections and of about 
10 6 events for the distributions. For the integrated cross sections, the numbers in the 
Tables quoted in parenthesis after the last digit correspond to the la MC errors. As for 
the distributions, each of them has been generated as an histogram consisting of 25 bins. 
Therefore, in the dominant region of a given distribution (e.g. in the central range of the 
W rapidity and lepton pseudorapidity and around the Jacobian peak of the W transverse 
mass and lepton transverse momentum), the statistical uncertainty associated to the cross 
section in a bin is at the level of some per mille. In the tails of the distributions, the MC 
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error can grow up to a few per cent level. This statistics is sufficient for a meaningful 
interpretation of the effects of interest in the present study and discussed in the following. 
Actually, whenever rather mild effects are registered (e.g. of a few per cent), they typically 
show up in the dominant regions where the statistical uncertainty is about a factor of ten 
smaller. On the other hand, in the tails of the distributions, characterized by a poorer 
statistics, the contributions of interest are typically quite large (at the level of several or 
tens of per cent) and are therefore meaningful in comparison with the size of the statistical 
fluctuations. Last but not least, the most moderate effects generally correspond, as it 
will be shown in the following, to the difference between the results of the additive and 
factorized recipes for the combination of EW and QCD corrections. However, since the 
two recipes are defined in terms of the very same theoretical ingredients, and use the same 
MC data, it turns out that predictions of the two combinations are strongly correlated 
variables and, as such, the uncertainty associated to each of these two predictions largely 
cancel out when considering the difference of the two. 

3.2 Monte Carlo tuning 

Before discussing the numerical results, we would like to spend some words about the effort 
done in MC tuning, in order to normalize all the codes to the same LO or NLO cross section 
and, therefore, avoid possible bias in the interpretation of the physical effects of interest. 

Concerning the program ResBos-A, we switched off some EW contributions present 
inside the code in the form of Improved Born Approximation and transformed the running 
W width in the W propagator into a fixed width, to normalize the code at a pure LO cross 
section. Moreover, we used the program without NLO QED final-state corrections when 
providing pure QCD predictions. 

To tune FEWZ against the other codes, we substituted the default choice of using the 
experimental branching ratio BR(PF — > [iv) with the tree-level theoretical branching ratio 
T(W — ► fj,v)/T(W — ► all) = 1/9, for consistency with the other used programs. 

As far as MC@NLO is concerned, we tuned it at the next-to-leading-order in QCD, by 
comparing its predictions with the corresponding NLO QCD results of FEWZ and MCFM. 
This tuning has been performed at the level of completely inclusive (i.e. without cuts) cross 
sections both for the Tevatron and the LHC, as well as in the presence of a cut on the 
invariant mass M\ v > 0.8 TeV at LHC energies, in order to check the correct normalization 
in the region important for new-physics studies. 



Monte Carlo 


ALPGEN 


FEWZ 


HORACE 


MCFM 


ResBos-A 


o"lo (pb) 


906.3(3) 


905.4(2) 


905.6(1) 


905.1(1) 


905.3(2) 



Table 3: MC tuning at the Tevatron for the LO cross sections of the sum of the processes pp — > 
W + — ► fi + + v and pp — > W~ — > fj,~ + P, according to the cuts of Table 1 and using the CT6EQ6M 
PDFs with scale M Xv . 

The results are shown in Table || for the LO cross sections at the Tevatron, in Table |] 
for the LO cross sections at the LHC, in Table || for the tuning between FEWZ, MCFM 



- 12 - 



Monte Carlo 


ALPGEN 


FEWZ 


HORACE 


MCFM 


o"lo (pb) 


8310(2) 


8306(1) 


8308(1) 


8305(1) 



Table 4: MC tuning at the LHC for the LO cross sections of the sum of the processes pp — > W + — > 
/i + + v and pp — > W~ — ► fj,~~ + D, according to the cuts of set up a. in Table 2 and using the 
MRST2004QED PDFs with scale y/p'i w + Mf v 



Monte Carlo 


FEWZ 


MC@NLO 


MCFM 


< v a tr ° n (pb) 


2635.5(4) 


2639.1(5) 


2640(1) 


^(Pb) 


21058(3) 


21031(3) 


21008(2) 



Table 5: MC tuning at the Tevatron and LHC for MC@NLO and MCFM NLO inclusive cross 
sections of the sum of the processes pp — > W~ —* fi~D and pp — > W + — > using the 

PDF sets CT6EQ6M at the Tevatron with scale M lv and MRST2004QED at the LHC with scale 
\fp\_w + M /t> respectively. 



Monte Carlo 


MC@NLO 


MCFM 


o"NLo(fb) 


50.34(1) 


50.28(2) 



Table 6: MC tuning at the LHC for MC@NLO and MCFM NLO cross sections of the sum of the 
processes pp — ► W~ — > pTv and pp — ► W + — ► in the presence of a cut on the invariant mass 
Mi v > 0.8 TeV and no other lepton selection criteria. The PDF set is MRST2004QED with scale 

VpIw + M l- 



and MC@NLO NLO inclusive cross section both at the Tevatron and LHC, and in Table 
|6| for the MCFM and MC@NLO predictions at LHC energies when imposing the invariant 
mass cut M\ v > 0.8 TeV. As can be seen, there is a very good agreement between the 
predictions of the different programs, since all the relative differences are at the few per 
mille level (or better), i.e. much smaller than the size of the radiative corrections discussed 



in the following. As already remarked in 15], the tuning procedure is essential because it 



validates the interpretation of the various contributions as due to physical effects and not 
to a mismatch in the set up of the codes under consideration. 

3.3 Integrated cross sections (Tevatron and LHC) 

In Table and Table [8| we present the results for the integrated cross sections correspond- 



ing to the selection criteria at the Tevatron and LHC quoted in Section 34, as obtained 
according to the following codes and theoretical recipes 

1. Leading order, as obtained by MCFM code in terms of the LO matrix element sup- 



plemented by the PDF sets defined in Section 3.1 



2. ALPGEN So: the hadron-level LO cross section of the CC DY process is convo- 
luted with parton shower evolution through HERWIG Parton Shower. Hence, such a 
prediction can be considered as analogous to that of a pure parton shower generator; 
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3. ALPGEN Si: the exact LO matrix elements for up to one parton radiation are 
interfaced to HERWIG parton shower, according to the MLM prescription of matrix 
elements-PS matching; 

4. ALPGEN S2: the exact LO matrix elements for up to two parton radiation are 
interfaced to HERWIG parton shower, according to the MLM prescription of matrix 
elements-PS matching; 

5. MCFM, i.e. pure NLO QCD predictions (MCFM NLO QCD), including also fac- 
torization and renormalization scale variation by varying them by a factor 1/2 or 2 



around the central values specified in Section |3j] (see Table 



6. HORACE NLO EW: HORACE, in the presence of NLO EW corrections; 

7. HORACE NLO+HERWIG PS: the NLO EW predictions of HORACE are interfaced 
with HERWIG QCD Shower evolution; 

8. MC@NLO: NLO QCD combined with HERWIG parton shower; 

9. ResBos-A CSS: ResBos-A, with pure QCD CSS resummation, without fixed-order 
perturbative QCD contributions and final-state QED radiation; 

10. ResBos-A NLO QED: without QCD resummation, without fixed-order perturbative 
QCD results but with final-state QED radiation; 

11. ResBos-A, as far as its complete predictions, including QCD resummation and final- 
state QED radiation, but without perturbative QCD corrections, are concerned; 



12. Equation ( |2.1| ), obtained by summing the predictions of MC@NLO with those of 
HORACE NLO convoluted with HERWIG PS according to eq. fl2l|) (additive com- 
bination) ; 



13. Equation fl2.2|) , for the combination of QCD and EW contributions in factorized form. 

As can be seen from Table [7| and, in particular from the comparison between the LO 
results and the ALPGEN So predictions, the QCD PS effects lower the cross section by 
about 10-r- 15% at the Tevatron and LHC set up a, while they are negligible at LHC set up 
b. The exact matrix element corrections, as present in ALPGEN S2, are not relevant for 
the cross section at the Tevatron, while they induce a some per cent (positive) contribution 
for the two set up at the LHC. 2 The NLO QCD corrections are positive, of about 8%, at 
the Tevatron, negative of about 2% at the LHC set up a and positive of about 20% at 
the LHC set up b. The EW corrections are at the level of a few per cent and negative 
at the Tevatron and LHC set up a., while they are quite large, of the order of —25%, for 
set up b. at the LHC, in agreement with the results of various studies in the literature 
on the virtual EW Sudakov logarithms. It is worth noticing that for the considered event 



2 The results of ALPGEN Si are not shown in Table |^ because they are practically indistinguishable, 
within the statistical errors of the MC integration, from the predictions of ALPGEN Sa. 
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Tevatron 


LHC 


Prediction 




a. b. 


Leading order 

ALPGEN S 

ALPGEN S 2 

MCFM NLO QCD 

HORACE NLO EW 

HORACE NLO+HERWIG PS 

MC@NLO 

ResBos-A CSS 

ResBos-A NLO QED 

ResBos-A 

Eq. Q QCD+EW add. 

Eq. (p|) QCD ® EW fact. LO 

Eq. (2.2) QCD ® EW fact. NLO 


905.1(1) 
815(1) 
814(1) 
979.1(6) 
881.5(2) 
792(1) 
967(3) 
933.4(4) 
877.9(3) 
920(3) 
944(3) 
925(3) 
915(3) 


8305(1) 7.128(3) 
7224(5) 7.15(1) 
7578(5) 7.29(1) 
8135(2) 8.681(2) 
8091(1) 5.569(2) 
7008(5) 5.55(1) 
8254(5) 8.64(1) 

8038(9) 7.04(2) 
7877(8) 6.71(2) 
7895(8) 6.50(1) 



Table 7: Integrated cross sections, relative to the selection criteria of Table 1 and Table 2, according 
to different approximations for the various programs considered in the study. For the Tevatron and 
set up a. at the LHC the cross sections are given in pb, while the units of measure are fb for set 
up b. at the LHC. 



selection the NLO QCD and EW corrections tend to compensate at the Tevatron and 
LHC set up b, while they sum up coherently at the LHC set up a. The QCD PS, when 
convoluted with the NLO EW contributions, decreases the cross sections of about 10%, as 
already noted above at the level of pure QCD results, at the Tevatron and the LHC set up 
a, while it is irrelevant at the LHC set up b. The combination of NLO QCD corrections 
with PS resummation, as predicted by MC@NLO, gives a correction at the per cent level 
w.r.t. the pure NLO QCD correction obtained with MCFM. The soft-gluon resummation 
a la CSS is responsible of a few per cent increase of the LO cross section at the Tevatron, 
while NLO final-state QED radiation is of the same order but of opposite sign, so that 
the two effects tend to compensate, yielding a full ResBos-A prediction of only about 1% 
larger than the LO cross section. On the other hand, the combination of EW and QCD 
corrections according to the (additive) formula of eq. (|2.1|) predicts an enhancement of 
the LO cross section at the Tevatron of about 7%, differing of some per cent from the 
ResBos-A result, as a consequence of the different QCD and EW ingredients between eq. 
( |2.1| ) and the ResBos-A formulation. The calculation of the combined EW and QCD as in 
the (factorized) formula of eq. ( |2.2j ) differs from the additive prediction of about 2%, thus 
reducing the difference w.r.t. ResBos-A. At the LHC, the combined QCD/EW correction 
is negative, of the order of few per cent, both in the set up a and b. It is worth noticing, 
however, that in the set up a this is due to the sum of mild corrections of the same sign, 
whereas in the set up b it is the result of the almost complete cancellation of large QCD and 
EW corrections. Also in this case, the difference of additive and factorized prescriptions 
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amounts to few per cent. 



Collider 


>/2) 
°NLO 


°NLO 


°NLO 


Tevatron 


978.5(8) 


979.1(6) 


987.1(6) 


LHC a. 


7875(1) 


8135(2) 


8393(2) 


LHC b. 


8.890(5) 


8.681(2) 


8.501(4) 



Table 8: NLO cross sections, including scale variation, as obtained by means of MCFM at the 
Tevatron and LHC, set up a. and b. The cross sections are given in pb for the Tevatron and set 
up a. at the LHC, in fb for set up b. at the LHC. 



In Table ^ we show the NLO predictions of MCFM at three different factorization/re- 
normalization scales, = = /V2,^,2/U, as it is customary, with /i = M\ v at the 



Tevatron and /i = y ^iv a ^ ^ ne LHC. The scale dependence at the Tevatron is at 

the 1% level, indicating a good stability of the NLO prediction for the assumed scales. At 
the LHC set up a., the cross section varies of about 3-4% as a function of the scale variation, 
while a 2% variation is observed for the cross section at the LHC set up b. As a whole, the 
results shown in Table [?] point out the following aspects: first of all, a proper combination 
of EW and QCD corrections is mandatory to achieve a cross section accuracy of the order of 
a few per cent, in particular for the set up b at the LHC where a subtle cancellation occurs; 
secondly, at the per cent accuracy the NNLO QCD corrections should be taken into account, 
being their contribution of a few per cent as estimated by the scale variations shown in 
Table |8| and confirmed by dedicated calculations in the literature 41, |l6|, |l7], [42]| ; last, 



at such a precision level also exact mixed 0(aa s ) contributions, at present not known, can 
be expected to be relevant as pointed out by the relative difference between the additive 
and factorized prescriptions of eqs. (24.) and (|2.2j), respectively. 



3.4 PDF uncertainties (LHC) 

In this Section we discuss how the PDFs uncertainties of experimental origin, as known of 
today, affect all the observables of interest for the physics analysis of the process pp — > fi+X 
at the LHC. 

In Figure || we show the results for the W rapidity and muon pseudorapidity (upper 
plots) and for the W transverse mass and muon transverse momentum (lower plots) accord- 
ing to set up a. specified in Table ^, as obtained with the NLO CTEQ61 parameterization 
available in the LHAPDF package. 

For each observable, we show the predictions corresponding to the maximum and 
minimum values returned by CTEQ61 PDFs, together with the result of the best fit parton 
densities. 

The PDF uncertainty for an observable O is the maximal change in O as a function of 
variables {z{} varying within the tolerance hypersphere |p3| : 

n/2 



i=l 
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CTEQfi 



CTEQ6 





Figure 2: CTEQ61 PDFs uncertainties for W rapidity and muon pseudorapidity (upper plots), 
and for the W transverse mass and muon transverse momentum (lower plots), according to set up a. 
at the LHC. In the lower panel of each plot, the relative deviations of the minimum and maximum 
predicted values w.r.t. the best fit PDF are shown. 



where 

xn. = T— ~ T 0{^i+t)-0{4-t) 
dzi ~ 2t 

n is the number of the error sets (30 for MRST, 40 for CTEQ), T is the tolerance, t is a 
small step in the space of z%. Here 0{z\, ± t, 2° / 2 ) is denoted as 0(zf ± t). 

In the lower panel of each plot, the relative deviations of the minimum and maximum 
predicted values w.r.t. the best fit PDF are shown. As can be seen, the spread of the 
predictions is at the level of some per cents for all the distributions. In particular, for the 
W rapidity and muon pseudorapidity the uncertainties vary within 2-4%, reaching their 
maximum value in the central rapidity region, while for the muon p± and the relative 
deviations are about 3-4% around the jacobian peaks, reaching the 6% level in the hard 
tails of the distributions. 
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Figure 3: The same as Figure 2 for MRST200fE parameterization. 

The same analysis is shown in Figure | for the NLO MRST2001E parameterization. 
The observed uncertainties are smaller by about a factor of two, as expected from the 



discussion addressed in Section |2T| about the different values of the tolerance parameter 
adopted in the two PDF sets. 

The results for the W transverse momentum are shown in Figure [|, both for CTEQ61 
(left plot) and MRST2001E (right plot). The uncertainties are almost flat in both cases, 
and of about 3% and 1.5% for CTEQ61 and MRST2001E parameterization, respectively. 

We also analyzed the PDF uncertainties on the transverse mass distribution when 
considering set up b. of Table || at the LHC, as shown in Figure || We observed that both 
the CTEQ61 and MRST2001E predictions display deviations w.r.t to the best fit at some 
per cent level from 1 to 2 TeV, while the uncertainties are, in the average, of the order of 
20% in the vicinity of 3 TeV, as due to the well-known gluon PDF uncertainties for large 
x values dominant in such a region. 

It is important to remind, as already remarked, that the estimate of PDF uncertainties 
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CTEQ6 



MRST2001 




Figure 4: CTEQ61 (left plot) and MRST2001E (right plot) PDFs uncertainties for the W trans- 
verse momentum distribution, according to set up a. at the LHC. 
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Figure 5: CTEQ61 (left plot) and MRST2001E (right plot) PDFs uncertainties for the W trans- 
verse mass in the high tail, according to set up b. at the LHC. 



obtained according to such a procedure are of experimental origin only, leaving aside other 
sources of uncertainty due to theory. 

At present the only available set which includes electromagnetic effects in the evolution 
of the PDFs is MRST2004QED. We have chosen this set to consistently subtract the 
initial state QED collinear divergences, and also to investigate the contribution of the 
extra photon- induced partonic subprocesses, due to the presence of a photon density in the 
proton. 



The recent PDFs sets [ 71 , |73|] have included new data and a refined treatment of the 
quark masses effects. These updates turned out to be important for DY phenomenology 
at the LHC: in fact the total cross sections within cuts change by approximately 3-6% 
w.r.t. the predictions obtained with the older sets. On top of the overall normalization 
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variation we observed a moderate (of the order of 1%) but not constant shape change when 
comparing the transverse mass or the lepton transverse momentum distributions computed 
with or without massive quark corrections. 

Our choice of using MRST2004QED, although it is not the most up-to-date proton 
parameterization, has been taken to keep the possibility of discussing the relevance of the 
photon-induced processes, which might be not negligible especially in the searches for new 
physics signals. For consistency at the Tevatron we have used the Resbos grids based on 
the CTEQ6.1 PDFs set, which does not include massive charm effects. 

The analysis of the present Section is not spoiled by the use of older sets of PDFs nor 
it does overestimate the uncertainties due to the experimental error from which the PDFs 
are affected. In fact we numerically checked that the 1-a spread of the predictions obtained 
with the most recent sets is of the same size, w.r.t. the results shown here. 

3.5 Numerical results for the Tevatron 

In this section, we present the results obtained for the process pp — ► W ± -> n± + X at 
the energy of the Tevatron Run II (y/s = 1.96 TeV), when imposing the cuts quoted in 
Table |. For the various observables of interest, we begin with a discussion of QCD effects, 
to continue with the analysis of the combination of EW and QCD corrections. The latter is 
realized according both to the additive formula of eq. Q2.1| ) and the factorized combination 
as in eq. fl2,2| ). The numerical results presented in the following have been obtained using 
HORACE including exact NLO EW corrections, but neglecting the contribution of higher- 
order effects due to multiple photon emission, since the phenomenological impact of these 
contributions has been already studied in detail, both at the Tevatron and LHC, in a series 



of previous papers [25, g^, 3C, 34]. Moreover, results for the W transverse momentum 
will not be shown because of the strong sensitivity of this observable to non-perturbative 
effects, that require a careful modeling and comparison with real data and are, therefore, 
beyond the scope of the present analysis. 

3.5.1 Observables for luminosity monitoring and PDF constraint 

After the successful MC tuning at the level of integrated cross sections discussed in Section 



3.2| , we performed detailed comparisons between the predictions of QCD codes under con- 
sideration for those observables of experimental interest to monitor the collider luminosity 
and to constrain the PDFs. 

In Figure |6| we show the results obtained for the W rapidity (left plot) and muon 
pseudorapidity (right plot) with ALPGEN S , ALPGEN S 2 , MC@NLO and ResBos-A 
CSS. 3 It is important to emphasize that in Figure || and in the following plots referring 
to distributions in the presence of QCD only, the results of all QCD programs have been 
normalized to the corresponding integrated cross section, in order to point out just the 
differences in the shape description. In the lower panels of Figure [6| we show the relative 
deviations of each QCD tool w.r.t. ResBos-A, chosen as a benchmark because of its wide 



3 As already stated for the integrated cross sections in Section 3.3, the results of ALPGEN Si are not 
shown in Figure ^ and in the next plots for the Tevatron because they are practically indistinguishable from 
the predictions of ALPGEN S2, as we explicitly checked. 
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vw >n 

Figure 6: W rapidity (left plot) and muon pseudorapidity (right plot) distributions according to 
the QCD predictions of ALPGEN S , ALPGEN S 2 , MC@NLO and ResBos-A. In the lower panels 
the relative deviations of each code w.r.t. ResBos-A are shown. 



use at the Tevatron 4 . It can be seen that, in spite of the different theoretical ingredients, 
the predictions of the QCD programs agree at the ~ 1% level almost in the whole shape. 
This can be easily understood because yw and r\\ are rather smooth distributions and, as 
such, quite insensitive to QCD shape differences, at least for the ranges accessible at the 
Tevatron. 

We also investigated the level of agreement between the QCD codes when considering 
their predictions for the W charge asymmetry, which is an important quantity at the 
Tevatron to derive information about the partonic contents of the proton. As can be seen 
from Figure |?], the absolute differences between the various predictions reach at most the 
1% level, again as a consequence of the quite smooth behavior of this distribution. 

A first illustration of the combination of EW and QCD corrections at the Tevatron is 
shown in Figure ^ for the yw and r\i distributions. The upper panels show the absolute 
predictions obtained by means of the codes ALPGEN So, MC@NLO and according to 
eq. ( |2.lD , when using MC@NL0 5 for the simulation of QCD effects in association with 
HORACE convoluted with HERWIG Parton Shower. In the lower panels, the relative 
effects in units of ALPGEN S* due to QCD (MCONLO), EW (HORACErerwig) and the 
combination of electroweak and QCD corrections (eq. fl2.1|) ) are shown, as obtained by 
appropriate combinations of the absolute predictions shown in the upper panel. Strictly 
speaking, the results referring to MC@NLO+HORACEherwig and denoted in the plots as 
EW correspond to the second term of eq. (2.1), in units of the predictions of ALPGEN Sq. 



It can be seen that EW corrections amount to a few per cent negative contribution, QCD 
NLO effects are of the order of 20%, resulting in a total combined correction of about 15%. 



4 More precisely, the code used at the Tevatron for the modeling of QCD in DY processes is ResBos. 
However, ResBos-A coincides with ResBos as far as QCD contributions are concerned. 

5 Instead of adopting MC@NLO, the same results could be obtained using ALPGEN with an overall 
jf-factor (defined as o~nlo /vlo), since the shape predictions of ALPGEN and MC@NLO are in good 
agreement, as previously shown. 
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Figure 7: The W charge asymmetry according to the QCD predictions of ALPGEN So, ALP- 
GEN Sa, MC@NLO and ResBos-A. In the lower panel the absolute deviations of each code w.r.t. 
ResBos-A are shown. 




Figure 8: W rapidity (left plot) and muon pseudorapidity (right plot) distributions according to 
MC@NLO, ALPGEN S and the additive EW+QCD combination. In the lower panels the relative 
effects due to QCD, EW and EW+QCD corrections are shown in units of ALPGEN S . 



We also studied the EW and QCD corrections to the W charge asymmetry and observed 
that they are of comparable size and tend to sum up, yielding a total correction at the 1% 
level on the average. We also performed a comparison with ResBos-A, to notice that the 
predictions of our EW/QCD recipe are in agreement with those of ResBos-A within about 
1% in the whole rji range. 

Figure |9] shows the comparison between the additive formula of eq. ( |2.1| ) and the factor- 
ized prescriptions of (|2.2|), as normalized to the LO hadron level distributions convoluted 
with HERWIG PS. The predictions lie in a few per cent band, the factorized prescriptions 
being quite close to one another. The relative difference between the additive and fac- 
torized recipes can be seen as an estimate of the uncertainty due to 0{a 2 s ) and 0{aa s ) 
corrections. 
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Figure 9: W rapidity (left plot) and muon pseudorapidity (right plot) distributions according to 
the additive EW+QCD combination and the two EW<g>QCD factorized prescriptions. In the lower 
panels the relative effects due to the different combinations are shown in units of ALPGEN Sq. 
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Figure 10: W rapidity (left plot) and muon pseudorapidity (right plot) distributions according to 
the additive EW+QCD combination, the two EW<g>QCD factorized prescriptions and ResBos-A. In 
the lower panels the relative differences of the various combinations w.r.t. ResBos-A are shown. 



In figure [ji] the absolute predictions of eqs. fl2.1| ) and ( |2.2| ) are compared with those 
of ResBos-A in the upper panel. As can be seen in the lower panels, the relative difference 
does not exceed the few per cent level, the agreement being better for the factorized 
prescriptions, consistently with the fact that ResBos-A adopts a factorized formulation 
as well. 

3.5.2 Observables for W precision physics 

The status of the QCD predictions for the observables relevant for precision measurements 
of the VF-boson mass is presented in Figure [ll|, showing the W transverse mass (left plot) 
and the muon transverse momentum (right plot), according to the results of ALPGEN So, 
ALPGEN S 2 , MC@NLO and ResBos-A, as in Figure || and Figure |7[ For such strongly 
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Figure 11: The same as Figure 6 for the W transverse mass (left plot) and muon transverse 
momentum (right plot) distributions. 



varying distributions, it can be noticed that the predictions of the QCD programs differ at 
some per cent level around the jacobian peak, which is the crucial region for M\y extraction. 
The relative differences can reach the 5-10% level in the distributions tails. Presumably, the 
discrepancies around 50 GeV for and 25 GeV for p l ± have to be ascribed to soft-gluon 
resummation, which is accurately described in ResBos-A but not taken into account at the 
same precision level in the other codes. On the other hand, the differences observed around 
90 GeV for and 50 GeV of p l j_ are probably due to hard collinear PS effects, which are 
incorporated in ALPGEN and in MC@NLO but are absent in ResBos-A. Actually, one can 
notice that in this region the predictions of ALPGEN S2 and MC@NLO well agree with 
the pure PS approximation of ALPGEN So, supporting the above interpretation. 




Mf (GeV) p' ± (GeV) 

Figure 12: W transverse mass (left plot) and muon transverse momentum (right plot) distributions 
according to MC@NLO, ALPGEN So and the additive EW+QCD combination. In the lower panels 
the relative effects due to QCD, EW and EW+QCD corrections are shown in units of ALPGEN Sq. 



A further example of the combination of EW and QCD corrections at the Tevatron 
is shown in Figure [l^ for the Mj^ and p l ± distributions. Similarly to Figure ^, the upper 



panels show the absolute predictions of ALPGEN So, MC@NLO and according to eq. (|2,1|), 
while the lower panels illustrate the relative effects due to QCD and EW corrections only 
and to the combination of electroweak and strong contributions. As can be observed, in 
the most interesting region around the jacobian peak the NLO corrections are partially 
compensated by the EW contributions, yielding a total correction of about 10-15%. It is 
worth noticing that the convolution of EW corrections with the QCD shower changes the 
shape of EW corrections themselves, smearing them w.r.t. a pure parton-level calculation. 




M ± (GeV) p' ± (GeV) 

Figure 13: W transverse mass (left plot) and muon transverse momentum (right plot) distributions 
according to the additive EW+QCD combination and the two EW<g>QCD factorized prescriptions. 
In the lower panels the relative effects due to the different combinations are shown in units of 
ALPGEN S . 



Figure 13 shows the comparison between the additive formula of eq. (2.1) and the 
factorized prescriptions of (|2.2j ), as normalized to the LO hadron level distributions con- 
voluted with HERWIG PS. Some comments are in order here. With respect to the case, 
already discussed, of yw and r/i, the differences between the two factorized prescriptions 
are sizable. In particular, in the case of p l j_, the factorized formula with LO normalization 
leads to pathological results. This can be understood as due to the fact that NLO correc- 
tions are dominated, in the hard p± tail, by 2 — > 3 QCD subprocesses not included in a 
LO calculation. For this kind of observables, the NLO calculation is the first perturbative 
order able to cover all the relevant phase space regions; hence it is the "lowest order" 
approximation w.r.t. which evaluate the size of radiative corrections. 

Figure 14 shows the comparison between the additive and factorized formulae on the 
one hand and the ResBos-A prediction on the other one. From the comparison with Fig- 
ure |ll] referring to QCD shape differences only, it can be seen that the trend of the relative 
deviations is very similar to that observed between the results of MC@NLO and ResBos-A, 
while the size of the differences increases of some per cents in the tails, as a consequence of 
the inclusion of the normalization in the codes predictions. It can also be noticed that the 
differences between the two approaches are largely dominated by QCD effects. Therefore, 
as already remarked about Figure [0], the differences below the peaks are probably due 
to soft-gluon resummation effects. On the other hand, the discrepancies present above 
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Figure 14: W transverse mass (left plot) and muon transverse momentum (right plot) distributions 
according to the additive EW+QCD combination, the two EW<g>QCD factorized prescriptions and 
ResBos-A. In the lower panels the relative differences of the various combinations w.r.t. ResBos-A 
are shown. 



the peaks have to be presumably ascribed to hard collinear PS contributions taken into 
account in our calculation but not implemented in ResBos-A. To better understand the 
origin of these differences, it would be interesting to perform comparisons with the predic- 
tions of ResBos-A when taking into account the effect of finite order QCD contributions 
(the so-called Y terms) in its formulation. Unfortunately, no public grids are available 
for such theoretical configuration. Concerning the EW contributions, the two programs 
essentially predict the same quantitative effects. This is not unexpected because EW cor- 
rections in the vicinity of the W peak are largely dominated by final-state QED radiation, 
as widely discussed in the literature, and therefore the inclusion of exact NLO corrections 
in HORACE versus the approximation limited to the inclusion of NLO final-state QED cor- 
rections in ResBos-A translates into differences at the 0.1% level, as we explicitly checked. 
As a whole, it is worth stressing that closely around the jacobian peak the additive recipe 
and ResBos-A show differences of a few per cent, for both and p l j_, whereas the factor- 
ized NLO prescription exhibits smaller deviations in the case. Such differences could 
be relevant in view of the future measurements of the W mass with ~ 20 MeV precision 
at the Tevatron and should be further scrutinized on the experimental side. In particular, 
the combination of MC@NLO with HORACE, as adopted in this study, could be used to 
cross-check the precision measurement of the H^-boson mass as currently obtained in terms 



of the generators ResBos (for QCD effects) and HORACE, PHOTOS @ and WGRAD 
(for EW corrections). 

Results for photon-induced processes at the Tevatron are not shown since they are 
negligible because of the very small photon content inside the proton for the typical x 
values probed by DY kinematics at the Tevatron. 

3.6 Numerical results for the LHC 

In this section, we present the results obtained for the process pp — » — > ^ + X at 
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the LHC (y/s = 14 TeV), when imposing the cuts quoted in Table As for the Tevatron 
analysis, we begin with a discussion of QCD effects on the relevant observables (considering 
shape predictions only), to continue with the study of the combination of EW and QCD 
corrections. Results taking into account higher-order leading log QED corrections and 
predictions for the W transverse momentum observable will not be given in the following, 



for the same reasons already emphasized at the beginning of Section 3.5 



3.6.1 Observables for luminosity monitoring and PDF constraint 



0.3 
0.25 

0.2 
0.15 

0.1 
0.05 

0.1 

0.05 





(i 

yw 



Figure 15: W rapidity (left plot) and muon pseudorapidity (right plot) distributions according 
to the QCD predictions of ALPGEN S , ALPGEN Si, ALPGEN S 2 and MC@NLO. In the lower 
panels the relative deviations of each code w.r.t. MC@NLO are shown. 

In Figure |l5| we show the results obtained for the W rapidity (left plot) and muon pseu- 
dorapidity (right plot) with the QCD generators ALPGEN S , Si and S 2 , and MC@NLO, 
considering set up a. in Table In the lower panels we show the relative deviations of 
the ALPGEN variants w.r.t. MC@NLO. It can be noticed that, while for r\\ the relative 
differences between ALPGEN and M C@NLO are at a few per cent level in the whole range, 
for the W rapidity the agreement in the shape predicted by the two generators is quite sat- 
isfactory (at the 1% level) in the central region, but it deteriorates in the very forward and 
very backward limits, that anyway marginally contribute to the integrated cross section. In 
particular, in these regions, the pure PS predictions represented by ALPGEN So indicate 
a quite smooth excess of events w.r.t. MC@NLO, while the (practically indistinguishable) 
deviations of ALPGEN Si and of ALPGEN S 2 w.r.t. MC@NLO have a more pronounced 
opposite trend. 



The interplay between EW and QCD corrections at the LHC is shown in Figure 16 for 
the yw and r\i distributions. The upper panels show the absolute predictions of MC@NLO, 
ALPGEN S , MCONLO+HORACE (additive combination) and MC@NLO®HORACE 
(factorized combination) interfaced to HERWIG PS. In the lower panels, the relative effects 
due to the combination of EW and QCD corrections, as well as for the EW and QCD 
contributions separately, are shown. The size of NLO EW and QCD corrections is almost 
the same as the one seen for the Tevatron in Figure |8[ In particular, the negative NLO EW 
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Figure 16: W rapidity (left plot) and muon pseudorapidity (right plot) distributions according to 
the predictions of MC@NLO, ALPGEN S , MC@NLO+HORACE and MC@NLO<X>HORACE. In 
the lower panels the relative effects due to QCD, EW and combined corrections are shown in units 
of ALPGEN S . 

effects partially cancel the positive contribution due to NLO QCD corrections, resulting 
into combined EW and QCD corrections at the 10% level, for both yyy and r\i and for 
both additive and factorized prescriptions. Hence, the interplay between EW and QCD 
contributions is crucial for precise simulations of the observables relevant for luminosity 
monitoring and PDF constraint at the LHC. 
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Figure 17: W rapidity distribution according to the predictions of MCFM, MC@NLO, 
MCFM+HORACE and MC@NLO+HORACE H erwig- In the lower panel the relative effects due 
to 0{a 2 s ) QCD corrections, 0(aa s ) mixed corrections and QCD+EW corrections with and without 
QCD showering. 

As remarked in Sect. 2.4, both the additive and factorized combination of EW and 
QCD corrections contain, beyond the very same 0(a) and 0(a s ) content, higher-order 
contributions dominated by 0(a 2 s ) and mixed 0(aa s ) corrections. It is therefore interesting 
to analyze how these higher-order corrections compare to each other. Since such higher- 
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order effects originate, in both the additive and factorized recipe, from the QCD PS and 
the convolution of the QCD PS with the NLO predictions, a simple strategy to disentangle 
the relative effect of 0(ot 2 s ) and 0(aa s ) corrections is to compare the predictions of the 
combination of pure NLO codes (lacking, by construction, the contribution of the QCD 
shower evolution) with those of one of the two recipes described in Sect. 2.4 or their 
variants. An example of such a study is given for yyy m Fig. [H], where, for defmiteness, 
the results of the additive recipe MC@NLO + HORACEherwig are compared with those 
of the pure NLO combination MCFM + HORACE, as well as with the predictions of the 
"intermediate" theoretical formulations given by MC@NLO + HORACE (missing mixed 
0(aa s ) corrections) and MCFM + HORACEherwig (missing 0{a 2 s ) shower effects). As it 
can be seen from the lower panel, the fixed NLO predictions of MCFM + HORACE differ 
from those of MC@NLO + HORACEherwig at the some per cent level, showing the overall 
importance of QCD shower and mixed EW-QCD corrections for a precise simulation of the 
yw distribution. It can be also noticed that 0(aa s ) corrections are very small, while QCD 
shower corrections amount to a few per cent, thus dominating the EW-QCD combination 
for this observable. 

3.6.2 Observables for W precision physics 






Figure 18: W transverse mass (left plot) and muon transverse momentum (right plot) distributions 
according to the QCD predictions of ALPGEN S , ALPGEN Si, ALPGEN S 2 and MC@NLO. In 
the lower panels the relative deviations of each code w.r.t. MC@NLO are shown. 

The QCD predictions for the W transverse mass and muon transverse momentum ac- 
cording to set up a. of Table || is shown in Figure [l8|. The results of different variants of 
the ALPGEN generator are compared with the predictions of MC@NLO. As can be seen 
from the lower panels, the agreement in the shape predicted by the two generators is very 
satisfactory (at the 1% level) around the jacobian peak, while the relative deviations reach 
the 5-10% level in the hard tail of both and p\ spectrum. In particular, it can be 
noticed that for the relative difference of ALPGEN S w.r.t. MC@NLO has the same 
size and shape of the deviations observed for the predictions of ALPGEN Si and of ALP- 
GEN S2, while this is not the case for the hard p l j_ tail, where the MC@NLO predictions 
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are in better agreement with the ALPGEN versions including matched matrix-elements 
corrections rather than with ALPGEN So- This can be easily understood in terms of the 
rather smooth dependence of the distribution from QCD corrections, at variance of 
the pronounced sensitivity of p l ± from hard QCD radiation. The importance of the latter 
is particularly visible in Figure 19, showing the predictions of the QCD generators for the 
muon transverse momentum distribution up to 80 GeV. Actually, one can see the well- 
known fact ]5l| that in the high tail of such distribution, i.e. above ~ 50 GeV, there is an 
important enhancement due to NLO and matched matrix elements corrections w.r.t. the 
QCD PS approximation represented by ALPGEN So- Furthermore, a substantial agree- 
ment between the shape predicted by MC@NLO and the ones obtained with ALPGEN Si 
and ALPGEN S2 is observed, although a closer inspection reveals discrepancies above 
~ 50 GeV in the 10% range, as illustrated in the lower panel of Figure [11], where however 
the cross section is quite small and the relative statistics quite limited. As a whole, the 
deviations between MC@NLO and ALPGEN in its different flavors shown in Figure 15 and 
Figure [l^ can be ascribed to the different matching procedures between real and virtual 
QCD radiation implemented in the two generators. 




Figure 19: The muon transverse momentum distribution up to 80 GeV according to the QCD 
predictions of ALPGEN S , ALPGEN Si, ALPGEN S 2 and MC@NLO. In the lower panels the 
relative deviations of ALPGEN Si and ALPGEN S 2 w.r.t. MC@NLO are shown. 



Because of the importance of precise predictions for the M\y determination from fits to 
Mj^ and p l ± , we show in Figure [20| the W transverse mass (left plot) and muon transverse 
momentum (right plot) distributions according to the QCD predictions of MC@NLO, when 
varying the renormalization/factorization scale from its default value /io = = = 



P±w + to /io/2 and 2[iq. As it can be seen from the lower panels in Figure 20, the 
scale variations induce relative differences w.r.t. the default choice at a few per cent level 
for both the distributions, without modifying the shape of the distributions themselves. 
These deviations reflect the ~ ±3% variations obtained for the cross sections by integrating 
the above distributions by MC@NLO, given by 8026(5) pb for 11 = fj, /2, 8254(5) pb for 
fi = fio and 8502(5) pb for fi = 2/iq. 
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Figure 20: W transverse mass (left plot) and muon transverse momentum (right plot) distributions 
according to the QCD predictions of MC@NLO for renormalization/factorization scale fi at its 
default value fio given in the text and for fj, = 1/2/zo, 2fiQ. In the lower panels the relative deviations 
due to scale variations w.r.t. the default choice are shown. 




Figure 21: The same as Figure 16 for the W transverse mass distribution (left plot) and muon 
transverse momentum (right plot). 



The combined effect of EW and QCD contribution at the LHC is illustrated in Figure 
21. The upper panels show the predictions of the generators ALPGEN So, MC@NLO, 
MC@NLO + HORACE and MC@NLO®HORACE. The lower panels illustrate the relative 
effects of NLO QCD and EW corrections, as well as their combination. From fig. [2l] it can 
be seen that the NLO QCD corrections are positive and tend to compensate the effect due 
to EW corrections. Therefore, their interplay is unavoidable for a precise Mjy extraction 
at the LHC, yielding an overall correction of about 5-10% close to the peaks for both the 
additive and factorized prescriptions. As already stressed for such distributions at the 
Tevatron, it is also worth noticing that the relative size and shape of EW corrections to 
M± and p l ± is modified by the convolution with QCD shower evolution, as can be inferred 
from comparison with the results for pure EW contributions existing in the literature. 
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Figure 22: Electroweak and 7-induced corrections without (left) and with (right) QCD shower 
evolution to the lepton transverse momentum distribution. 

This indicates that a naive combination of QCD and EW corrections, in the absence of 
QCD shower evolution on top of EW effects, is rather inadequate for precise simulation of 
physical observables. This feature is remarked in Figure showing the relative effect of 
EW corrections and photon-induced processes for the lepton p± distribution in the absence 
of QCD parton shower convolution (left plot) and in the presence of QCD shower evolution 
(right plot). In particular, it can be noticed that the some per cent enhancement due to 
the 7-induced processes in the hard tail of the p l j_ distribution without QCD shower largely 
disappears when the EW effects are convoluted with QCD shower radiation. 
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Figure 23: The same as Fig. 17 for the W transverse mass distribution. 
An inspection of the relative impact due to 0(a^) and mixed 0(aa s ) corrections analo 



gous to that shown in Fig. |17] for yyy is shown in Fig. |23| for the transverse mass distribution. 
For this observable, it can be seen that the QCD shower evolution is particularly important 
below the Jacobian peak, where it introduces corrections of some per cent, while mixed 
EW-QCD contributions are not negligible above it. More importantly, a naive combina- 
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tion of QCD and EW corrections, omitting the contribution of QCD shower, appears to 
be inadequate for a precision calculation of Mjf , as the clearly visible difference between 
the predictions of MC@NLO + HORACE HE rwig and of MCFM + HORACE points out. 
It follows immediately that a measurement of the W mass at the LHC with an aimed ac- 
curacy of a few MeV must necessarily rely upon MC generators including combined QCD 
and EW corrections according to the highest theoretical standards. 6 

3.6.3 Observables for new-physics searches 

The numerical results corresponding to the distributions interesting for the search for new 



physics at the LHC are shown in Figure 24 and Figure 25 for what concerns QCD predic 



tions, and in Figure 26 and Figure E7], for the combination of EW and QCD corrections. 



The simulations have been performed using the cuts referring to set up b. of Table |2] and 
imposing additional jet veto conditions for what concerns the results shown in Figure |2?]. 

Concerning Figure ^4], it can be seen that the shape predicted by the different variants 
of ALPGEN w.r.t. the one obtained with MC@NLO has a very similar trend for both 
Mj^ in the range 1-2 TeV and p l ± in the region 0.5-1 TeV. The predictions of all the 
ALPGEN versions, both without and with matrix element corrections, nicely agree with 
the MC@NLO results within a 5% band. 





Figure 24: The same as Figure 18 for the transverse mass (left plot) and lepton transverse mo- 
mentum (right plot) distributions in the very high tails, according to set up b. at the LHC. 



Figure ^5] shows the W transverse mass (left plot) and muon transverse momentum 
(right plot) distributions according to the QCD predictions of MC@NLO, when varying the 
renormalization/factorization scale from its default value /Uo = fJ*R = I^f = ■J P±w + 
to Ho/2 and 2/xrj. As can be seen from the lower panels, the scale variations induce relative 
differences w.r.t. the default choice of a few per cent around 1 TeV for , reaching 

6 It is worth noticing that for the first bins close to the kinematical boundary [M^ — 50 GeV, not 
shown in the figure) the prediction by MCFM is negative. This is a well known effect already discussed 
in p2| , Q, and due to perturbative instabilities of the NLO calculation. As a consequence, the results by 
MCFM around the jacobian peak appear to be slightly larger than the corresponding ones by MCCDNLO, 
without contradicting the results for the integrated cross sections shown in Tab. M. 
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O(-10%) at 3 TeV, but where the cross section is two orders of magnitude smaller. For p]_ 
the scale variations are milder, always within ±5%. 




Figure 25: The same as Figure 20 for the transverse mass (left plot) and lepton transverse mo- 
mentum (right plot) distributions in the very high tails, according to set up b. at the LHC. 



The interplay between QCD and EW corrections in the high tail of and p l ± 

distributions is shown in Figure ^6| (for the standard cuts of set up b. in Table |2|) and 
Figure p7| (when including additional jet veto conditions). 
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Figure 26: The same as Figure 21 for the transverse mass (left plot) and lepton transverse mo- 
mentum (right plot) distributions in the high tails, according to set up b. at the LHC. 



The jet veto conditions have been implemented by vetoing any hard jet with p± > 
30 GeV in the central rapidity region (|j/jet| < 2.5). In both cases, and for both 
and p l ± , NLO QCD corrections are positive (of the order of 20-40% for standard cuts and 
between 20-60% in the presence of jet veto requirements) and combine with very large 
negative EW Sudakov logarithms. When imposing standard selection criteria, their fac- 
torized combination is about -10(+10)% for Mf ot 1.5(3) TeV and about -5(-20)% for 
p 1 ^ ~ 0.5(1) TeV, the additive combination yielding smaller effects. Therefore, a precise 
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normalization of the SM background to new physics searches necessarily requires the si- 



multaneous control of QCD and EW corrections, as also pointed out in [15| for the NC DY 
process in the high invariant mass region, where a similar compensation between strong 
and EW contributions takes place. It is worth noticing that in such a region the large 
negative Sudakov virtual logarithms can be compensated by real weak boson emission, as 
pointed out in [95]. Concerning the inclusion of two-loop EW Sudakov logarithms, whose 
calculation is available in the literature [96, 97, 98 1, it should be considered in view of a 
theoretical accuracy at the percent level. For the same reason it should be accompanied 
by the simultaneous control of the NNLO QCD corrections. 
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Figure 27: The same as Figure 21, according to set up b. at the LHC and in the presence of 
additional jet veto conditions. 




Figure 28: The same as Figure 17 for the transverse mass (left plot) and lepton transverse mo- 
mentum (right plot) distributions in the high tails, according to set up b. at the LHC. 



As done for y\y and MjT in the presence of the standard analysis cuts denoted as 
LHCa, we also investigated how the predictions of a combination of pure NLO codes for 
EW and QCD corrections compare with those in which QCD shower effects are taken into 
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account when considering Mjf and p l T in the set up LHC b. The results of this study are 
shown in Fig. 25 for Mjf (left plot) and p l T (right plot). It can be seen that the pure NLO 
results differ from those including QCD shower contributions at the 10% level for pL and 
even more, between 20-30%, for Mjf in the very hard tail. This important difference is 
due to the presence of large EW Sudakov logarithms in the NLO EW calculation, giving 
rise to sizable 0(aa s ) corrections at the NNLO level. Hence, the evaluation of the DY 
background to new physics searches in this region should take carefully under control the 
combined effect of QCD-EW corrections, especially in the event new physics would manifest 
as a moderate modification of the DY continuum, rather than as a clearly visible resonance. 
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Figure 29: Upper panels: transverse mass distribution of the single W production process (black 
histogram) and of the ti process (red histogram), according to the event selection specified in the 
text. Lower panels: relative difference of the two contributions, for two transverse mass windows. 



To complete the phenomenological analysis, we also performed an investigation of the 
contribution to the transverse mass due to the top-pair production process pp — > ti — > 
^ + X (simulated with ALP GEN) m comparison with pp -> W ± -> /z± + X. The 
obtained results are shown in Figure [2^, for a moderate transverse mass range up to 700 
GeV (left plot) and in the high tail up to 3 TeV (right plot). For the ti process, all 
the t leptonic decays are simulated and the events are considered as contributions to the 
signature when at least a muon is present in the final state. Whenever just a single muon 
is part of the final state products, the W transverse mass is reconstructed in terms of its 
transverse momentum, whereas in the case two muons are present the W transverse mass 
is calculated in terms of the leading-j^ muon. The unobserved X state consists of the two 
6-quarks and two neutrinos coming from the top and W decays, respectively. The missing 
transverse momentum entering the definition of W transverse mass is obtained by summing 
over the missing pt of the two neutrinos and requiring for each of the two the condition 
$T > 25 GeV, consistently with the cuts imposed in the analysis of W production events. 
As can be seen from Figure 2£, the ti process gives, for the event selection assumed, a 
sizable (of the order of several tens per cent) contribution to the p ± + X signature in the 
transverse mass range between about 100 GeV and 600 GeV, while it is at a few per cent 
level in the high tail. 
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3.7 Results for jet multiplicity (Tevatron and LHC) 



We conclude our phenomenological study with the presentation of the results concerning 
the distributions of the number of jets, showing the predictions obtained for the Tevatron 
in comparison with those valid for the LHC. The jets are required to satisfy the following 

< 5, ARij > 0.7, where Ai?j ? - is the separation between the 



plane. Jets are reconstructed by means the cone algorithm 



cuts: E Tj > 20 GeV, |%| 
i-th and j-th jet in the r/ - 
provided by the GET JET package [98], which represents a simplified jet cone algorithm a 
la UAl (for the sake of simplicity, we stop the evolution at the shower level). 
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Figure 30: QCD predictions of ALPGEN S , ALPGEN Si, ALPGEN S 2 and MC@NLO for the 
number of jets distribution at the Tevatron (left plot) and at the LHC, set up a. (right plot). In 
the lower panels the relative deviations of each code w.r.t. MC@NLO are shown. 



From the QCD generators used in our analysis, the predictions for the number of jets 
distribution are shown in Figure for the Tevatron (left plot) and the LHC (right plot), 
set up a. At the Tevatron energies, about 90% of single W production events are without 
any extra real hard QCD radiation, while the fraction of events with one additional jet is 
about 10%. At the LHC, the sharing of the events is, not unexpectedly, quite different: 
about 80% of the events do not contain extra jets, the fraction of events with one additional 
jet is slightly smaller than 20%, while events with al least two additional jets amount to a 
few per cent. As it can be seen from the lower panels in Figure |30|, there is good agreement 
between the predictions of MC@NLO and ALPGEN variants for the W production rate 
without any extra jet both at the Tevatron and LHC. At the Tevatron, for the events 
with one additional jet, MC@NLO predictions differ from ALPGEN So results (where the 
additional jet is generated through the PS cascade) at the 5% level and from ALPGEN Si 
and ALPGEN S 2 predictions at the ~ 10% level. At the LHC, all the ALPGEN variants 
provide results differing from the MC@NLO ones at the ~ 10% level for events with one 
extra jet, while for the fraction of two additional jet events there is good agreement between 
MC@NLO and ALPGEN Si but significant deviations are present between MC@NLO and 
either ALPGEN S or ALPGEN S 2 , as expected. 

In Figure 31 we show the number of jets distribution corresponding to set up b. at the 



-37- 



1 

0.8 

I 0.6 

Jit 0.4 

0.2 


0.4 

0.2 



-0.2 
-0.4 



MCSNLQ - 

ALPGEN S a - 

ALP GEN Si ■ 

ALPGEN Sa " 



Figure 31: QCD predictions of ALPGEN S , ALPGEN Si, ALPGEN S 2 and MC@NLO for the 
number of jets distribution at the LHC, set up b. In the lower panels the relative deviations of each 
code w.r.t. MC@NLO are shown. 



LHC. As can be seen, the requirement MW > 1 TeV changes very significantly the number 



of jets distribution w.r.t. the situation of the more inclusive set up a. shown in fig. 30 



In particular, the fraction of events without extra jets is about 30% and it's exceeded by 
the fraction of events with one additional jet, at the 40% level, at least for what concerns 
the ALPGEN Si and ALPGEN S2 predictions. There are important fractions of events 
containing two and three extra jets, at the 20% and 10% level, respectively. From the lower 
panel, it can be noticed that the relative differences between the results of ALPGEN Si 
and ALPGEN S2 and those of MC@NLO lie in a ±20% range for the predictions relative 
to the fraction of events with at least one extra jet. 



4. Conclusions and perspectives 

In this paper we presented a detailed phenomenological analysis of the EW and QCD 
contributions, as well as of their combination, to single W boson production in hadronic 
collisions, considering the energies of interest for the experiments at the Fermilab Tevatron 
and the CERN LHC and in the presence of realistic selection criteria. In our study, we 
made use of MC tools which include the recent advances in EW and QCD calculations. 

We analyzed all the quantities of interest for the many facets of W physics programme 
at hadron colliders and showed that a far from trivial interplay between EW and QCD 
effects is present for most of the observables. We noticed, in particular, partial cancel- 
lations between EW and QCD corrections for the experimental distributions of interest 
for both precision studies (such as the luminosity monitoring or the measurement of the 
W mass) and searches for new physics. Remarkably, in the region of the high transverse 
mass tail above 1 TeV, important for the search for new gauge bosons, these cancellation 
occurs almost completely between huge positive QCD corrections and very large negative 
contributions due to EW Sudakov logarithms. This emphasizes the need for a careful com- 
bination of strong and EW contributions in present and future analyses of the CC DY 
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process, as also pointed out in [15] for the NC DY channel in the high invariant mass 
region at the LHC. We also remarked that the convolution of the EW effects with QCD 
shower evolution is definitely relevant for a correct simulation of the distributions, since 
the relative size and shape of EW contributions is considerably modified when compared 
with the same features in the absence of the combination with a QCD PS. This can be 
understood in terms of the modifications introduced by QCD PS on the kinematics of the 
final-state leptons w.r.t. a pure EW calculation. We also discussed in detail and evaluated 
the uncertainty inherent the combination of EW and QCD corrections by comparing dif- 
ferent theoretical recipes (additive vs factorized) and showing that they can differ at a few 
per cent level. 

In relation to the accuracy of the theoretical tools presently available and used by 
the experimental collaborations in comparison with the precision already reached at the 
Tevatron and foreseen at the LHC, some final comments are in order. While the predictions 
of the EW programs are in very good agreement and the situation can be considered well 
under control (also when taking into account multiple photon emission), the same can 
not be said, strictly speaking, for the QCD generators currently used in the experimental 
analyses at the Tevatron or for simulations by ATLAS and CMS collaborations. Actually, 
we noticed, for example, differences of 2% or more in predictions for the shape of the 
distributions relevant for W mass extraction between ResBos-A and our combination of 
EW and QCD tools. Therefore, it is our opinion that a cross-check of the Tevatron W mass 
results, generally obtained by means of ResBos in association with WGRAD or PHOTOS, 
with the predictions of our recipe, based on the combination of MC@NLO with HORACE, 
would be highly desirable, in view of a robust, high-precision measurement of M\y, which 
is such an important input for indirect constraints on the Higgs boson mass. On the other 
hand, we concluded that the results of the different QCD generators are in very good 
agreement for what concern the W rapidity and lepton pseudorapidity distributions, and 
that EW corrections are well under control for such distributions. This reinforces the need 
for deeper attempts to monitor the hadron collider luminosity with presently available tools 
in terms of such distributions, and this goal should be attainable at the Tevatron with a 
few per cent precision, a certainly useful "exercise" to pave the way to LHC collaborations 
for future luminosity measurements along this direction. 

Concerning the forthcoming data taking at the LHC, we agree with the conclusions 



of (Tq, 17| about the feasibility of using the process-independent module PHOTOS to 



simulate photonic corrections around the W peak for early stages of data analysis, but we 
also provided clear evidences that such an approach to the treatment of EW corrections will 
not be sufficient in the later stages of analysis at high integrated luminosity. This caveat 
applies to high-precision extraction of the W boson properties at the LHC, as well as to 
predict correctly the SM background to new-physics searches in the high tail, where pure 
EW Sudakov logarithms and not photonic effects dominate within the full set of one-loop 
corrections. In particular, for such a region, two-loop Sudakov contributions and NNLO 
QCD corrections would be needed, together with the contribution of real W, Z emission. 

In general, for the LHC we remarked that available calculations and tools do not 
currently allow to reach a theoretical accuracy better than some per cent level, when 
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excluding PDF uncertainties. If this could be acceptable for earlier stages of analysis, future 
measurements at the LHC would probably require the calculation of complete 0{aa s ) 
corrections. All the necessary theoretical ingredients are not for the time being implemented 
into a single generator, which would represent the optimal solution for simulation and 
analysis of DY process in hadronic collisions, in its many and quite different aspects. 

For readers' convenience, we summarize in a final Table the relative effects of the 
different sources of corrections to the integrated cross section. We also include a further 
Table reporting an estimate of the theoretical accuracy of the best predictions that combine 
QCD and EW corrections according to the different recipes proposed. In tab. || NLO QCD 





NLO QCD 


NLL QCD 


NLO EW 


Shower QCD 


0(aa s ) 


Tevatron 


8 


16.8 


-2.6 


-1.3 


~ 0.5 


LHC a 


-2 


12.4 


-2.6 


1.4 


~ 0.5 


LHC b 


21.8 


20.9 


-21.9 


-0.6 


~ 5 



Table 9: Relative effect of the main sources of QCD, EW and mixed radiative corrections to the 
integrated cross sections for the Tevatron, LHC a and LHC b. 



is the complete 0(a s ) correction, NLL QCD is the matrix element contribution of the 
NLO QCD correction, NLO EW is the full 0{a) correction, Shower QCD stands for the 
0(a'g),n > 2 correction and 0(aa s ) represents the mixed EW-QCD corrections estimated 
by properly combining the additive and factorized cross sections. It is worth noticing in 
particular that the latter corrections remain below the 1% level for typical event selections 
at the Tevatron and the LHC, while they can amount to some per cent in the region 
important for new physics searches at the LHC. 



*(%) 


5a /a (scale) 


5a /a (FA) 


5a 1 a 


Tevatron 


~ 1 


~ 2 


2 


LHC a 


~ 2.5 


~ 2 


2.5 


LHC b 


~ 1.5 


~ 5 


5 



Table 10: Estimate of the present theoretical accuracy for the calculation of the integrated cross 
section at the Tevatron, LHC a and LHC b. 



Concerning the theoretical accuracy of the integrated cross sections calculation shown 
in Tab. 1C, it is a measure of the missing/incomplete higher order o? s and aa s contributions. 
It has been assessed by neglecting the effect of the PDF uncertainties and according to the 
following procedure: the relative scale variation of the additive cross section 5a /a (scale) 
and the relative difference between the additive and factorized cross sections 5a /a (FA) 
have been computed and compared; as error estimate, the largest of the two entries has 
been taken. 

Some comments are in order here. The above error estimate are in agreement with phe- 
nomenological results available in the literature about the size of NNLO QCD corrections, 
that are known to contribute at the ~ 2% level for standard cuts at the Tevatron and the 
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LHC. The errors quoted in Tab. 1C do not include the contributions of the two- loop EW 
Sudakov logarithms, that amount to a few per cent limited to the very high W transverse 
mass/lepton transverse momentum tails. Last, the error estimate reported can be con- 
sidered as rather conservative. Actually, the factorized cross sections used in the present 
analysis contain the bulk of NNLO QCD contributions, as can be inferred by comparing the 
relative difference between the additive and factorized predictions with the exact results 
for NNLO corrections given in [38, 39, fE], 41] and finding them in pretty fine agreement. 
Hence, a more aggressive estimate of the theoretical error of the factorized formulae could 
be derived by taking the relative difference the factorized formulae, obtaining 0(1%) at 
the Tevatron and the LHC a and 0(3%) at LHC b. This estimate could be put on firmer 
grounds through detailed comparisons with the NNLO exact QCD calculations available 
in the literature. Furthermore, in the very hard tails of the distributions, the calculation 
of still unavailable aa s corrections should be performed. 

Possible perspectives of the present work would be a careful analysis of how the preci- 
sion measurement of the W boson mass and width is affected by the combination of EW 
and QCD corrections here proposed and an application of the same approach to the study 
of the NC DY channel. 



Acknowledgments 

We are indebted to all the colleagues who allowed us, with their invitations, to present 
preliminary results of our work at various conferences and workshops: M. Czakon, H. 
Czyz, G. Degrassi, S. Ferrag, S. Forte, J. Gluza, M. Grazzini, S. Jadach, W. Placzek, M. 
Srzypeck, D. Wackeroth, Z. Was and G. Zanderighi. We thank K. Ellis and J.M. Campbell 
for help with the use of the code MCFM and valuable correspondence, S. Frixione for helpful 
comments on the use of the code MC@NLO and M.L. Mangano for useful discussions and 
interest in our work. We acknowledge M.R. Whalley for help on LHAPDF. We are also 
grateful to the colleagues of the common paper |15[| for fruitful collaboration. We wish to 
thank several experimental colleagues, in particular M. Bellomo, G. Polesello, A. Tricoli 
and V. Vercesi, for useful discussions and interest in our work. Useful discussions with 
P. Nason during the INFN Workshop on Monte Carlo's, Physics and Simulations at the 
LHC []10Cf are gratefully acknowledged. 

The work of A. V. was supported by the European Community's Marie-Curie Research 
Training Network under contract MRTN-CT-2006-035505 (HEP-TOOLS). 



References 



1] D. Denegri, Phys. Rept. 403 (2004) 107 and references therein 



[2] V.M. Abazov et al., CDF and D0 collaborations, \Phys. Rev. D 70 (2004) 09200g , 



hep-ex/0311039; 



CDF Collaboration and D0 Collaboration, Combination of CDF and DO results on the W 
boson mass and width, arXiv : 0808 . 0147 [hep-ex]; 



T. Aaltonen et al., CDF Collaboration, [Phys. Rev. D 77 (2008) 112001] , arXiv : 0708 . 3642 
[hep-ex] . 



- 41 - 



[3] 
[4] 
[5] 

[6] 
[7] 



D. Acosta et al., \Phys. Rev. D 71 (2005) 052002| , friep-ex/041105S 



D. Acosta et al., \Phys. Rev. D 71 (2005) 051104 frep-ex/0501025 . 

U. Baur, Electroweak physics at the Tevatron and LHC: Theoretical status and perspectives, 



hep-ph/0511064 



D. Wackeroth, Electroweak Physics: Theoretical Overview, tiep-ph/0610058 



S. Haywood et al., Electroweak Physics, in Proceedings of the Workshop on Standard Model 
Physics (and more) at the LHC, G. Altarelli and M.L. Mangano eds., CERN Report 2000-04, 
p. 117. 



W.T. Giele and S. Keller, Phys. Rev. D 57 (1998) 4433, iep-ph/970441S 



V. Bulge, Ch. Jung, G. Quast, A. Ghezzi, M. Malberti, T. Tabarelli de Fatis, Prospects for 
the precision measurement of the W mass with the CMS detector at the LHC, CMS AN 
20006/033; 

N. Besson, M. Boonekamp, E. Klinkby, T. Petersen and S. Mehlhase, for ATLAS 



Collaboration, Eur. Phys. J. C 57 (2008) 627, arXiv : 0805 . 2093 [hep-ex] 



[10] 

[11] 
[12] 

[13] 

[14] 
[15] 

[16] 

[17] 

[18] 

[19] 
[20] 
[21] 
[22] 
[231 



D. Froidevaux and V.A. Mitsou, Experimental prospects at the Large Hadron Collider, 
arXiv: 0905. 0258 [hep-ex]. 



M. Dittmar, F. Pauss and D. Zurcher, \Phys. Rev. D 56 (1997) 7284 |hep-ex/9705004 



V.A. Khozc, A.D. Martin, R. Orava and M.G. Ryskin, \Eur. Phys. J. C 19 (2001) 313 



hep-ph/0010163 . 

W.T. Giele and S.A. Keller, Hard scattering based luminosity measurement at hadron 



collider, liep-ph/0104053 



G. Polesello and M. Prata, Eur. Phys. J. C 32S2 (2003) 51 



C. Buttar et al., Standard Model Handles and Candles Working Group: Tools and Jets 
Summary Report, Les Houches 2007, Physics at TeV colliders, arXiv : 0803 . 0678 [hep-ph] . 



N. E. Adam, V. Halyo, S. A. Yost and W. Zhu, |j. High Energy Phys. 0809 (2008) 133] , 
arXiv: 0808. 0758 [hep-ph]. 



N. E. Adam, V. Halyo and S. A. Yost, \J. High Energy Phys. 0805 (2008) 062j , 
arXiv: 0802. 3251 [hep-ph]. 

G. Balossini et al., J. Phys. Con}. Ser. 110 (2008) 042002; 

G. Balossini et al., [Acta Phys. Polon. B 39 (2008) 16751, arXiv: 0805. 1129 [hep-ph]: 



G. Balossini et al., \Acta Phys. Polon. B 38 (2007) 3407 



D. Wackeroth and W. Hollik, \Phys. Rev. D 55 (1997) 6788| , frep-ph/9606398 



U. Baur, S. Keller and D. Wackeroth, \Phys. Rev. D 59 (1999) 013002| , |hep-ph/9807417 



S. Dittmaier and M. Kramer, Phys. Rev. D 65 (2002) 073007, hep-ph/0109062 



U. Baur and D. Wackeroth, Phys. Rev. D 70 (2004) 073015, hep-ph/0405191 



V. A. Zykunov, \Eur. Phys. J. Direct. C 3 (2001) | , |hep-ph/0107059| ; 
V. A. Zykunov, Phys. Atom. Nucl. 71 (2008) 732; 
V. A. Zykunov, Phys. Atom. Nucl. 69 (2006) 1522. 



- 42 - 



A. Arbuzov, D. Bardin, S. Bondarenko, P. Christova, L. Kalinovskaya, G. Nanava and 
R. Sadykov, \Eur. Phys. J. C 46 (2006) 407| , |hep-ph/0506110| . 



CM. Carloni Calame, G. Montagna, O. Nicrosini and A. Vicini, J. High Energy Phys. 12 



(2006) OH, hep-ph/0609170. 



S. Brensing, S. Dittmaier, M. 1. Kramer and A. Muck, \Phys. Rev. D 77 (2008) 073006 , 
arXiv: 0710. 3309 [hep-ph] . 

C. Buttar et al., Les houches physics at TeV colliders 2005, standard model and Higgs 
working group: Summary report, hep-ph/0604120. 

C.E. Gerber et al., Tevatron-for-LHC Report: Top and Electroweak Physics, 
arXiv : 0705 . 325 1 [hep-ph] . 



C. M. Carloni Calame, G. Montagna, O. Nicrosini and M. Treccani, Phys. Rev. D 69 (2004) 



037301, hep-ph/0303102. 



C. M. Carloni Calame, G. Montagna, O. Nicrosini and M. Treccani, J. High Energy Phys. 
0505 (2006) 0191 ^iep-ph/050221£ 



W. Placzek and S. Jadach, Eur. Phys. J. C 29 (2003) 325, hep-ph/0302065. 



P. Golonka and Z. Was, \Eur. Phys. J. C 45 (2006) 97j |hep-ph/0506026 . 

K. Hamilton and P. Richardson, \j. High Energy Phys. 0607 (2006) 010| , |hep-ph/0603034 . 



CM. Carloni Calame, S. Jadach, G. Montagna, O. Nicrosini and W. Placzek, Acta Phys 



Polon. B 35 (2004) 1643| , |hep-ph/0402235 



G. Altarelli, R. K. Ellis and G. Martinelli, Nucl. Phys. B 157 (1979) 461 



R. Hamberg, W. L. van Neerven and T. Matsuura, Nucl. Phys. B yVucl. Phys. B 359 (1991) 
JH [Erratum ibid. 644 (2002) 403]. 



J. M. Campbell and R. K. Ellis, Phys. Rev. D 65, 113007 (2002), hep-ph/0202176 ; 
J. M. Campbell and R. K. Ellis, Phys. Rev. D 62, 114012 (2000), hep-ph/0006304 ; 
J. M. Campbell and R. K. Ellis, Phys. Rev. D 60, 113006 (1999), hep-ph/9905386. 



C. Anastasiou, L. J. Dixon, K. Mclnikov and F. Petriello, \Phys. Rev. Lett. 91 (2003) 182002| , 
hep-ph/0306192. 



C. Anastasiou, L. J. Dixon, K. Melnikov and F. Petriello, \Phys. Rev. D 69 (2004) 094008| , 
hep-ph/0312266. 



K. Melnikov and F. Petriello, Phys. Rev. Lett. 96 (2006) 231803, tiep-ph/0603182 



K. Melnikov and F. Petriello, Phys. Rev. D 74 (2006) 114017, tiep-ph/0609070 



S. Catani, L. Cieri, G. Ferrera, D. de Florian, M. Grazzini, Vector boson production at 
hadron colliders: a fully exclusive QCD calculation at NNLO, arXiv : 0903 . 2120 [hep-ph] . 



Phys. Rev. Lett. 103 (2009) 082001 



C. Balazs and C. P. Yuan, \Phys. Rev. D 56 (1997) 5558| , |hep-ph/9704258 



F. Landry, R. Bro ck, P.M. Nadolsky and C.-P. Yuan, ^hys. Rev. D 67 (2003) 073016| , 
hep-ph/0212159 . 



[45] S. Frixione and B. R. Webber, J. High Energy Phys. 0206 (2002) 029, tiep-ph/0204244. 



-43- 



P. Nason, \,J. High Energy Phys. 0411 (2004) 0401, [hep-ph/0409146|; 

S. Frixione, P. Nason and C. Oleari, \j. High Energy Phys. 0711 (2007) 070 , 
arXiv: 0709. 2092 [hep-ph] . 



S. Alioli, P. Nason, C. Oleari and E. Re, \J. High Energy Phys. 0807 (2008) 060 , 
arXiv: 0805. 4802 [hep-ph]. 



M.L. Mangano, M. Moretti, F. Piccinini, R. Pittau and A.D. Polosa, J. High Energy Phys 
0307 (2003) 001| , jiep-ph/0206295 . 



A. Schalicke and F. Krauss, \J. High Energy Phys. 0507 (2005) 01§ [hep-ph/0503281 ; 

T. Gleisberg et. al., \J. High Energy Phys. 0402 (2004) 056| , |hep-ph/0311263] ; 

T. Gleisberg et al., \j. High Energy Phys. 0902 (2009) 007| , arXiv: 0811 . 4622 [hep-ph]. 



T. Stelzer and W.F. Long, Comput. Phys. Commun. 81 (1994) 357, hep-ph/9401258; 



F. Maltoni and T. Stelzer, J. High Energy Phys. 02 (2003) 02, tiep-ph/0208156 



J. Alwall et. al., |j. High Energy Phys. 0709 (2007) 028| , arXiv: 0706. 2334 [hep-ph]. 



A. Kanaki and C. G. Papadopoulos, Comput. Phys. Commun. 132 (2000) 306; 



A. Cafarella, C. G. Papadopoulos and M. Worek, Comput. Phys. Commun. 180 (2009) 1941 



P.M. Nadolsky, AIP Conf.Proc. 753 (2005) 158, tiep-ph/0412146 



U. Baur, Int. J. Mod. Phys. E 17 (2008) 826, hep-ph/0701164. 



P.Z. Skands, AIP Conf. Proc. 792 (2005) 73, |hep-ph/0507129 . 



D.A. Kosower, AIP Conf. Proc. 842 (2006) 189, |hep-ph/060401E 



Q.-H. Cao and C.-P. Yuan, \Phys. Rev. Lett. 93 (2004) 04200l| , |hep-ph/0401026 



S. Jadach, M. Skrzypek, P. Stephens, Z. Was and W. Placzek, Acta Phys. Polon. B 38 



(2007) 2305. 



S. Frixione and M.L. Mangano, J. High Energy Phys. 05 (2004) 056, tiep-ph/0405130 



H. Jung et al., Proceedings of the workshop: HERA and the LHC workshop series on the 
implications of HERA for LHC physics, arXiv : 0903 . 3861 [hep-ph] ; 

S. Alekhin et al., HERA and the LHC: A Workshop on the implications of HERA for LHC 
physics: Proceedings Part A and Part B, |hep-ph/0601012| and |hep-ph/0601013| . 

R.S. Thorne, A.D. Martin, R.G. Roberts and W.J. Stirling, AIP Conf. Proc. 792 (2005) 365, 
hep-ph/0507015 . 

M.R. Whalley, D. Bourikov, R.C. Group, The Les Houches accord PDFs (LHAPDF) and 



LHAGLUE, iep-ph/0508110 



A.D. Martin, R.G. Roberts, W.J. Stirling and R.S. Thorne, \Eur. Phys. J. C 28 (2003) 455| , 
hep-ph/0211080. 



J. Pumplin, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky and W.K. Tung, J. High Energi 



Phys. 0207 (2002) 012, liep-ph/0201195 



D. Stump, J. Huston, J. Pumplin, W. K. Tung, H. L. Lai, S. Kuhlmann and J. F. Owens, \J.\ 
High Energy Phys. 0310 (2003) 046] , |hep-ph/0303013 . 



[64] A. Tricoli, A. Cooper-Sarkar and C. Gwenlan, Uncertainties on W and Z production at the 
LHC, |hep-ex/0509002 



-44- 



A.D. Martin, R.G . Roberts, W.J. Stirling and R.S. Thorne, \Eur. Phys. J. C 35 (2004) 325| , 
hep-ph/0308087 . 

J. Huston, J. Pumplin, D.R. Stump and W.K. Tung, \J. High Energy Phys. 0506 (2005) 080| , 
hep-ph/0502080. 



A. D. Martin, R. G. R oberts, W. J. Stirling and R. S. Thorne, \Eur. Phys. J. C 39 (2005) 



155} [hep-ph/0411040 



M. Roth and S. Weinzierl, \Phys. Lett. B 590 (2004) 190[ |hep-ph/0403200 



H. Spiesberger, \Phys. Rev. D 52 (1995) 4936J , |hep-ph/9412286 



J. Kripfganz and H. Perlt, Z. Physik C 41 (1988) 319 



A.D. Martin, W.J. Stirling, R.S. Thorne, G. Watt, \Phys. Lett. B 652 (2007) 292| , 
arXiv : 090 1 . 0002 [hep-ph] . 



P. M. Nadolsky et al., \Phys. Rev. D D78 (2008) 013004 [arXiv: 0802 . 0007 [hep-ph] . 
R.S. Thorne, A.D. Martin, W.J. Stirling, G. Watt, arXiv : 0907 . 2387 [hep-ph]. 



G. Corcella et al., HERWIG 6.5 release note, liep-ph/0210213 



T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, S. Mrenna and E. Norrbin, Comput 



Phys. Commun. 135 (2001) 23S, hep-ph/0010017. 



R. Kuhn, F. Krauss, B. Ivanyi and G. Soff, Comput. Phys. Commun. 134 (2001) 223 



hep-ph/0004270 



F. Krauss, A. Schalicke and G. Soff, Comput. Phys. Commun. 174 (2006) 876 



hep-ph/0503087 



S. Gieseke, P. Stephens and B. Webber, J. High Energy Phys. 12 (2003) 04E 



hep-ph/0310083 



T. Sjostrand and P. Skands, Eur. Phys. J. C 39 (2005) 129, hep-ph/0408302 



T. Sjostran, S. Mrenna and P. Skands, \Comput. Phys. Commun. 178 (2008) 852| , 
arXiv : 07 10 . 3820 [hep-ph] . 



F. Caravaglios and M. Moretti, \Phys. Lett. B 358 (1995) 332j, [hep-ph/9507237 ; 

F. Caravaglios, M.L. Mangano, M. Moretti and R. Pittau, \Nucl. Phys. B 539 (1999) 215| , 



hep-ph/9807570 



S. Catani, F. Krauss, R. Kuhn and B.R. Webber. \J. High Energy Phys. 11 (2001) 063 



hep-ph/0109231 



L. Lonnblad, \J. High Energy Phys. 05 (2002) 046| , |hep-ph/0112284 



M. Mangano, The so-called MLM prescription for ME/PS matching , talk given at the 
Fermilab ME/MC Tuning Workshop, October, 2004. 



S. Mrenna and P. Richardson, \J. High Energy Phys. 05 (2004) 040| , [hep-ph/0312274 



M.L. Mangano, M. Moretti, F. Piccinini and M. Treccani, \J. High Energy Phys. 01 (2007) 



013j , [hep-ph/0611129 



S. Hoche et al, Matching parton showers and matrix elements, pep-ph/0602031 



J. Alwall et al., Eur. Phys. J. C 53 (2008) 473, arXiv : 0706 . 2569 [hep-ph] 



J.C. Collins, D.E. Soper and G. Sterman, Nucl. Phys. B 250 (1985) 199 



-45 - 



K.-P. O. Diener, S. Dittmaier and W. Hollik, \Phys. Rev. D 72 (2005) 093003 
hep-ph/0509084 . 



W. Hollik, T. Kasprzik and B. A. Kniehl, \Nucl. Phys. B 790 (2008) 13g , arXiv: 0707. 2553 
[hep-ph] . 



J. H. Kuhn, A. Kulesza, S. Pozzorini and M. Schulze, \Nucl. Phys. B 797 (2008) 2% 
arXiv: 0708. 0476 [hep-ph] ; 



J. H. Kuhn, A. Kulesza, S. Pozzorini and M. Schulze, \Phys. Lett. B 651 (2007) 160 
hep-ph/0703283 . 

A. Denner, S. Dittmaier, T. Kasprzik and A. Muck, Electroweak corrections to W 
hadroproduction including leptonic VF-boson decays, arXiv : 0906 . 1656 [hep-ph] . 



jet 



E. Maina, S. Moretti and D. A. Ross, Phys. Lett. B 593 (2004) 143 [Erratum-ibid. B 614 



(2005) 216], hep-ph/0403050 



P.M. Nadolsky and Z. Sullivan, PDF uncertainties in WH production at Tevatron, 
hep-ph/0110378 . 

S. Catani and B. R. Webber, JHEP 9710 (1997) 005 [arXiv:hep-ph/9710333] and references 
therein. 



U. Baur, \Phys. Rev. D 75 (2007) 013005| , |hep-ph/0611241 



A. Denner, B. Jantzen and S. Pozzorini, \Nucl. Phys. B 761 (2007) l| , |hep-ph/0608326 . 



B. Jantzen, J. H. Kuhn, A. A. Pcnin and V. A. Smirnov, \Nucl. Phys. B 731 (2005) 188 
[Erratum-ibid. B 752 (2006) 327], [hep-ph/0509157. 



B. Jantzen, J. H. Kuhn, A. A. Pcnin and V. A. Smirnov, Phys. Rev. D 72 (2005) 051301 



[Erratum-ibid. D 74 (2006) 019901], bep-ph/0504111 



F.E. Paige and S.D. Protopopescu, in Physics of the SSC, Snowmass, 1986, Colorado, edited 
by R. Donaldson and J. Marx. 

[100] F. Ambroglini et al., Proceedings of the Workshop on Monte Carlo's, Physics and 
Simulations at the LHC PART I, arXiv : 0902 . 0293 [hep-ph] ; 
F. Ambroglini et al., Proceedings of the Workshop on Monte Carlo's, Physics and 
Simulations at the LHC PART II, arXiv : 0902 . 0180 [hep-ph] . 



- 46 - 



